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THE  DESIGN  OF  OBSERVERS 
FOR  THE  MATCALS  SYSTEM 

ABSTRACT 

Three  observers  are  designed  for  a  reduced  order  system  that  represents 
the  lateral  system  of  the  F4J  aircraft  in  an  automatic  landing  configura¬ 
tion.  The  observers  are  to  be  used  in  the  aircraft's  lateral  control  sys¬ 
tem  to  estimate  its  lateral  position  and  lateral  velocity,  in  place  of  the 
a-B  filter  that  is  currently  used  to  estimate  position  and  velocity.  Re¬ 
sults  that  are  obtained  from  simulations  of  the  F4J  aircraft  lateral  control 
system  indicate  that  an  observer  may  be  used  to  improve  the  system's  re¬ 
sponse. 
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I.  INTRODUCTION 


This  report  gives  the  results  of  the  continuation  of  the  studies  des- 
scribed  in  the  interim  report  of  this  contract  [1],  The  interim  report 
described  the  design  of  a  observer  to  replace  the  a-8  filter  in  the  flight 

dynamics  and  control  module  of  the  F4J  aircraft  lateral  control  system. 

This  system  is  a  part  of  the  Marine  Air  Traffic  Control  and  Landing  System 
(MATCALS). 

The  observer  systems  are  compared  to  four  different  a-8  systems.  The 
first  system  is  the  SPN-42  system  [3],  which  is  called  the  a-B  system  in 
this  report.  The  second  system  is  the  SPN-42  system  with  the  controller 
gains  reduced  by  50%  to  reduce  the  system  noise  response.  The  third  sys¬ 
tem  is  called  the  ydot  system,  and  is  described  in  Chapter  II.  The  fourth 

system  is  called  the  model  2  system,  and  has  been  used  in  mode  1  (automatic) 
landings  [3].  The  criteria  of  comparison  for  the  systems  are  the  radar- 
noise  response  and  the  wind  response. 
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II.  MATCALS  SYSTEM 


This  chapter  first  contains  a  brief  description  of  the  lateral  con¬ 
trol  system  of  the  F4J  aircraft  in  an  automatic  landing  configuration 
(MATCALS).  Next  various  system  design  criteria  are  developed.  Finally, 
two  variations  of  this  control  system  are  presented.  For  the  first  vari¬ 
ation,  the  loop  gain  is  reduced  by  50  percent.  In  the  second  variation, 
it  is  assumed  that  the  aircraft  lateral  velocity  is  available  as  a  feed¬ 
back  signal . 


F4J  Lateral  Control  System 

A  complete  description  of  the  F4J  aircraft  in  the  MATCALS  configura¬ 
tion  is  given  in  [2].  A  brief  description  of  this  system  will  be  given 
here. 

A  block  diagram  of  the  F4J  lateral  control  system  is  given  in  Figure 
2-1.  The  block  labeled  "aircraft  system"  is  the  aircraft  with  autopi¬ 
lots.  The  signal  $(t)  is  the  command  input  to  the  bank  autopilot,  and 
y(t)  is  the  lateral  distance  of  the  aircraft  from  the  extended  centerline 
of  the  runway.  The  signal  w{t)  models  the  wind  input  to  the  aircraft. 

The  signal  "radar  noise"  represents  the  measurement  inaccuracies  in  the 
determination  of  the  aircraft  position  by  the  radar,  and  yR(k)  is  the 
aircraft  position  as  measured  by  the  radar.  The  a-B  filter-controller 
combination  is  basically  a  proportional-integral-derivative  (PID) 
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Aircraft  Lateral  Control  System  with  the  «-$  Filter 


controller  used  to  compensate  the  closed-loop  system.  This  combination 
also  contains  filtering  to  reduce  the  effects  of  the  radar  noise  signal. 


System  Design  Criteria 

Even  though  the  MATCALS  control  systems  are  discrete  in  nature,  the 
important  design  criteria  can  be  developed  and  understood  by  considering 
the  equivalent  analog  system.  The  equivalent  analog  system  is  considered 
here  because  the  derivations  are  much  simpler. 

The  analog  equivalent  of  the  lateral  control  system  of  Figure  2-1 
is  given  in  Figure  2-2.  In  the  controller,  Kj  is  the  gain  in  the  inte¬ 
gral  path,  Kp  is  the  gain  in  the  proportional  path,  and  KQ  is  the  gain 
in  the  derivative  path.  For  the  F4J  lateral  control  system,  Kp  *  Kp  - 
0.75,  and  Kj  =  .00333  at  close  range.  The  noise-filtering  transfer  func¬ 
tion  has  been  omitted,  since  it  does  not  affect  the  following  derivations. 
The  transfer  function  from  bank  command  input  to  position  output  for  the 
aircraft  is  G-j(s),  and  from  wind  input  to  position  output  is  G2 ( s } .  The 
aircraft  velocity  y(t)  is  needed  for  later  derivations.  Yrn(s)  ‘•s  the 
radar,  noise  signal,  and  Ycom(s)  is  the  system  command  signal. 

Two  transfer  functions  of  critical  importance  to  system  operation 
are  the  transfer  function  from  the  radar  noise  signal  to  the  position 
output  and  the  one  from  the  wind  to  the  position  output.  From  Figure 
2-2,  we  see  that 

YfO  (Kp+KT/s+Kns)G1(s) 

Y^TiT  "  T T  (yKj/s+KpSjfi^'s)  (2-]) 

and 
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It  is  necessary  that  the  plane  response  very  little  to  the  radar 


noise.  However,  the  transfer  function  for  the  command  input,  Ycorn>  to 
position  output  is,  from  Figure  2-2, 


V  (si 

Y 

com 


(Kp+iys+KpS^Cs) 

TsT  =  i  +  (Kp+iys+^slG^s) 


(2-3) 


This  transfer  function  is  identical  to  that  of  Y(s)/Y  (s).  Thus,  if 
system  response  to  the  radar  noise  is  reduced,  the  system  response  to 
the  system  command  is  also  reduced.  For  low  frequencies; 


(Kp+KI/s+KDs)G1(s) 


»  1 

s=jw 


(2-4) 


Then,  from  (2-1)  and  (2-3), 


m. 

Ycom 


(sT 


S=j<d 


(2-5) 


Hence,  to  reduce  noise  response,  we  must  effectively  reduce  system  band¬ 
width;  i.e.,  we  must  reduce  the  frequency  range  over  which  (2-4)  is 
satisfied. 

However,  reducing  the  frequency  range  over  which  (2-4)  is  satisfied 
increases  the  system  response  to  wind.  If  (2-4)  is  satisfied,  then  for 
this  frequency  range  (2-2)  becomes 


This  gain  is  small,  since  the  denominator  is  large.  But,  if  the  denomi¬ 
nator  is  decreased  in  order  to  reduce  the  response  to  radar  noise,  the 
response  to  wind  is  increased. 

In  summary,  the  lateral  control  system  must  be  redesigned  in  such  a 
manner  that  the  system  response  to  radar  noise  is  reduced,  while  the 
open-loop  gain  as  given  in  (2-4)  remains  large  such  that  system  wind  re¬ 
sponse  is  not  degraded. 

System  with  Reduced  Gain 

The  first  attempt  to  satisfy  the  above  design  specification  is  sim¬ 
ply  to  reduce  the  controller  gains,  Kp,  Kj,  and  KQ  in  Figure  2-2,  by  50 
percent.  This  is  obviously  not  a  complete  solution,  since  the  system 
wind  response  is  increased.  It  has  been  suggested  that  this  gain  reduc¬ 
tion  will  give  an  acceptable  radar  noise  response,  and  for  this  reason 
this  system  will  be  used  here  for  comparison  purposes.  However,  the  wind 
response  for  this  case  will  probably  be  excessive,  and  thus  is  not  accept 
able.  It  is  then  assumed  that  this  system  does  give  an  acceptable  radar 
noise  response,  and  an  unacceptable  wind  response.  The  study  of  the  sys¬ 
tem  with  full  gains,  i.e.,  the  SPN-42  system,  gives  an  unacceptable  radar 
noise  response  and  an  acceptable  wind  response. 

The  system  with  full  gains  will  be  referred  to  as  the  a- 6  system, 
and  the  system  with  the  gains  reduced  by  50  percent  will  be  referred  to 
as  the  half-gain  system.  Thus  the  redesign  problem  is  to  design  a  system 
that  approaches  the  noise  response  of  the  half-gain  system  while  approach 
ing  the  wind  response  of  the  a-S  system.  Of  course,  adequate  stability 
margins  must  also  be  maintained  in  the  redesign. 
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A  second  variation  on  the  a-B  system  that  is  useful  to  study,  but 
cannot  be  implemented  because  of  hardware  constraints,  is  the  system  of 


Figure  2-3.  This  system  will  be  referred  to  as  the  ydot  system.  In 
this  system  it  is  assumed  that  the  aircraft  velocity  y(t)  is  measured 
directly  on  the  aircraft,  and  is  then  transmitted  to  the  ground-based 
controller.  Thus  the  differentiation  of  a  noisy  signal  in  the  PID  con¬ 
troller  is  no  longer  required,  and  this  leads  to  a  reduced  noise  response. 

For  the  ydot  system  of  Figure  2-3,  it  is  seen  that  the  following 
transfer  functions  apply. 


Y(s)  =  Y(s)  (KP+KI/s)G1(s) 

Y(s)  _  G2^s) 

Wti)  "  1  +  (Kp+KI/s+KDs)G1 (s) 


(2-7) 

(2-8) 


A  comparison  of  (2-7)  for  this  system  to  (2-1)  for  the  a-B  system  indi¬ 
cates  that  the  bandwidth  of  the  system  from  the  radar  noise  input  to 
position  output  has  been  reduced.  Furthermore,  comparing  (2-0)  to  (2-2) 
for  the  a-B  system  shows  that  the  wind  response  is  the  same  for  the  two 
systems.  Thus  the  ydot  system  appears  to  achieve  the  design  objectives, 
and  as  will  be  shown  later,  does  in  fact  achieve  the  design  objectives. 


f 
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III.  OBSERVER  CONTROL  SYSTEM 


A  possible  improvement  in  the  MATCALS  closed-loop  system  is  to  re¬ 
place  the  o-e  filter  in  Figure  2-1  with  a  different  filter.  The.  purpose 
of  the  a-6  filter  is  to  produce  a  filtered  aircraft  position  and  velocity 
signal,  y(k)  and  y  respectively,  based  on  the  radar  data.  An  observer, 
or  state  estimator,  is  a  different  type  of  filter  that  can  be  used  to 
produce  estimates  of  the  states  of  the  aircraft.  Observer  theory  is 
covered  in  detail  in  [1],  and  will  not  be  presented  here.  The  observer- 
based  MATCALS  lateral  control  system  is  shown  in  Figure  3-1.  This  system 
will  be  referred  to  as  the  observer  system. 

Note  in  Figure  3-1  that  the  observer  has  as  inputs  the  radar  data 
yR(k)  and  the  controller  output  $(k).  The  observer  outputs  are  the  es¬ 
timated  aircraft  position  y(k)  [x^(k)  in  the  observer  equations  below] 
and  estimated  aircraft  velocity  [^(k)].  Thus  there  is  an  additional 
loop  in  the  observer  system  that  does  not  appear  in  the  o-6  system,  which 
should  lead  to  increased  flexibility.  The  observer  designed  for  the  F4J 
lateral  control  system  will  now  be  presented. 

Observer  Design 

The  observer  designed  for  the  F4J  lateral  control  system  is  based 
on  a  third-order  approximation  to  the  ninth-order  model  for  the  aircraft 
system  in  Figure  3-1  [1],  This  discrete  model  is  given  by 
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of  the  F4J  Aircraft  Lateral  Control  System  with  the  Observer. 


‘x^k+1)1 

*1.0  0.1  0.00477116* 

x-|  (k) 

*0.000114237* 

x2(k+l) 

0.0  1.0  0.0932144 

CM 

X 

+ 

0.00338736 

x3(k+l) 

0.0  0.0  0.867429 

__ 

lx3^. 

|_0. 0661790 

<f>(k) 


y(k)  =  [1.0  0.0  0.0] 


x.j  (k) 

x2(k) 

x^(k) 


(3-2) 


In  these  equations,  x^(k)  =  y(k)  [aircraft  position],  x2(k)  =  y(k)  [air¬ 
craft  velocity],  and  x3(k)  =  y(k)  [aircraft  acceleration].  These  equa¬ 
tions  are  of  the  general  form 


x(k+l )  =  Ax(k)  +  B<j>(k) 
y(k)  =  Cx(k) 

The  equations  of  the  observer  are  given  by 

x(k+l)  =  (A-LC)x(k)  +  B$(k)  +  Ly(k) 

where  the  matrices  A,  B,  and  C  are  given  in  (3-3),  (3-1),  and  (3-2).  The 
matrix  L  is  determined  by  the  design  procedure.  To  determine  L,  the 
characteristic  equation  of  the  observer,  denoted  by  o(z),  must  be  chosen. 
Then,  from  (3-4), 

a(z)  =  |zl  -  (A-LC)|  =  0  (3-5) 

The  only  unknown  in  (3-5)  is  L,  and  Ackermann's  Formula  [1]  may  be  used 
to  solve  (3-5)  for  L. 

For  the  observer  designed  in  [1],  a(z)  was  chosen  as 


(3-3) 


(3-4) 


c(z)  -  (z-0.8)3 


(3-6) 


which  resulted  in 
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0.467428 


(3-7) 

r  is  the  root  time  constant,  then  roots  at 
time  constant  for  the  observer  of 


In  order  to  determine  the  effects  of  the  choice  of  a(z)  on  the  lat¬ 
eral  control  system  of  Figure  3-1,  observers  were  also  designed  for 

a(z)  =  (z-0.85)3  (3-8) 

and 

a(z)  =  (z-0.75)3  (3-9) 

The  observer  time  constant  for  (3-8)  is  then  0.615  seconds,  and  for  (3-9), 
0.348  seconds.  The  results  of  this  design  are  tabulated  in  Table  3-1. 

A  comparison  of  the  responses  of  these  three  observer  systems  to  the 
three  systems  described  in  Chapter  II  will  be  given  in  the  next  chapter. 


L  *  |  0.57863 
I  0.0353145 


Since  z  =  esT  *  e"T^T,  where 
z  *  0.8  in  (3-6)  result  in  a 


-0.1/t  _  Q  Q 


or 


t  *  0.448  seconds 
since  T  is  0.1  seconds. 
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TABLE  3-1 
Observer  Ma trices 

(z-.85)3  (z-.8)3 


0.317429 

0.467428 

0.254144 

0.57863 

0.00060367 

0.0353145 

(z-.75)3 

0.617429 

1.04757 

0.186368 


IV.  COMPARISON  OF  SYSTEMS 

In  this  chapter  a  comparison  is  made  of  the  radar  noise  responses 
and  the  wind  responses  of  the  F4J  lateral  control  system  for  the  six  dif¬ 
ferent  control  configurations  described  in  Chapters  II  and  III.  The  con¬ 
figurations  are: 

1.  a-e  system. 

2.  half-gain  system. 

3.  ydot  system. 

4.  (z-0.85)  observer  system  (characteristic  equation). 

o 

5.  (z-0.8)  observer  system. 

3 

6.  (z-0.75)  observer  system. 

3 

Recall  that  the  (z-0.85)  observer  is  the  slowest  one  ( i . e . ,  narrowest 
bandwidth)  and  the  (z-0.75)  observer  is  the  fastest  one  (broadest  band¬ 
width). 


Initial  Condition  Response 

The  first  comparison  of  the  systems  is  between  the  initial  condition 
responses.  All  responses  discussed  in  this  report  were  obtained  from  the 
simulation  given  in  [1]. 

For  the  initial  condition  responses,  the  radar  noise  signal  and  the 
wind  Input  were  set  to  zero.  The  initial  aircraft,  displacement  y(0)  was 
set  to  twenty  feet  at  a  range  of  13,223  feet,  which  required  a  flight 
time  to  touchdown  of  sixty  seconds.  The  results  are  given  in  Figures 
4-1  and  4-2. 
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y(t) ,  feet 


Figure  4-1.  Initial  Condition  Response  of  MATCALS  Systems. 


-B  system 


Figure  4-2.  Initial  Condition  Response  of  MATCALS  Systems. 


In  Figure  4-1  only  the  ( Z-.85 )'3  observer  system  response  is  plotted; 

3 

the  responses  of  the  other  two  observer  systems  lie  between  the  (z- . 85 ) 
observer  system  response  and  the  a-B  system  response.  The  ydot  system 
response  is  much  slower  than  the  a-8  system  response,  since  the  ydot  sys¬ 
tem  bandwidth  is  reduced  as  shown  by  (2-7).  The  half-gain  system  response 
is  also  slower  than  the  a-B  system  response,  because  of  the  reduced  sys¬ 
tem  bandwidth. 

Conclusions:  the  observer  systems  do  not  degrade  the  system  time 
response  to  the  extent  of  those  of  the  ydot  and  half-gain  systems.  How¬ 
ever,  the  time  response  of  the  half-gain  system  may  be  adequate;  this  is 
not  known. 


Stability  Margins 

An  important  measurement  of  the  relative  stability  of  a  closed-loop 
system  is  the  stability  margins,  i.e.,  the  gain  margin  and  the  phase  mar¬ 
gin. 

The  gain  margins  and  the  phase  margins  of  the  six  systems  are  given 
in  Table  4-1.  Recall  from  Figure  3-1  that  the  observer  system  contains 
two  loops;  thus  the  stability  margins  must  be  calculated  for  the  system 
opened  in  each  loop.  A  property  of  observer  systems  is  that  stability 
margins  do  not  change  in  the  primary  loop,  which  in  Figure  3-1  is  the 
system  opened  at  $(k).  However,  stability  margins  may  degrade  when  opened 
at  the  input  to  the  plant,  which  in  Figure  3-1  is  the  signal  <f>(t). 

Note  in  Table  4-1  that  stability  margins  for  the  observer  systems 
opened  at  $(k)  are  different  from  those  of  the  a-B  system.  This  differ¬ 
ence  originates  in  the  use  of  a  third-order  model  of  the  aircraft  system 
to  design  the  observers. 
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TABLE  4-1 

Stability  Margins-Opened  at  <$(k) 


system 

phase  margin 

gain  margin 

o-B 

49° 

15  dB 

bal f-gai n 

42° 

21  dB 

ydot 

49° 

15  dB 

(z-0.85)3 

55° 

11.2  dB 

(z-0.8)3 

55° 

11 .7  dB 

(z-0.75)3 

55° 

12.0  dB 

Opened  at  <t>(t) 

(z-0.85)3 

43° 

10.5  dB 

(z-0.8)3 

45° 

12  dB 

(z-0.75)3 

48° 

13.4  dB 
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Conclusions:  it  appears  that  the  stability  margins  of  all  systems 
are  adequate,  since  the  half-gain  system  has  the  smallest  phase  margin. 
Note  that  the  stability  margins  improve  as  the  observe  speed-of-response 
increases. 


Closed-loop  Noise  Frequency  Response 
Shown  in  Figure  4-3  are  the  closed-loop  frequency  responses  Y/Yrn 
for  the  six  systems.  The  signal  Y  is  the  aircraft  position,  and  Yrn  is 
the  radar  noise  signal.  Note  that  the  bandwidth  of  the  a-B  system  is 
the  broadest,  and  the  ydot  system  has  the  narrowest  bandwidth. 

Conclusions:  The  bandwidths  of  the  observer  systems  are  somewhat 
smaller  than  that  of  the  a-B  system,  with  the  (z- . 85 )  observer  (the 
slowest  one)  having  the  narrowest  bandwidth.  Thus  this  observer  system 
should  have  the  least  response  to  radar  noise,  when  compared  to  the 
other  two  observer  systems  and  the  a-B  system.  The  bandwidth  of  the  ydot 
system  is  approximately  one- third  that  of  the  a-B  system.  The  bandwidth 
of  the  half-gain  system  is  approximately  one-half  that  of  the  a-B  system. 

Closed-loop  Wind  Frequency  Response 
Given  in  Figure  4-4  is  the  closed-loop  frequency  response  Y/W,  where 
Y  is  the  aircraft  position  and  W  is  the  wind  input.  The  frequency  re¬ 
sponse  for  the  ydot  system  is  approximately  the  same  as  that  of  the  a-B 
system  and  could  not  be  plotted  as  a  separate  curve.  This  result  is  ex¬ 
pected,  as  explained  in  Chapter  2. 

Conclusions:  The  a-B  system  frequency  response  is  acceptable,  and 
that  of  the  half-gain  system  is  not.  The  observer  systems  frequency 
response  falls  between  these  two,  with  the  (z-.75)  observer  system 


4-6 


4-8 


I 

(the  fastest)  closest  to  the  acceptable  a-B  system  response  Tr.is  result 
is  expected,  since  the  fastest  observer  system  has  the  broadest  open-loop 
bandwidth,  and  this  is  the  requirement  for  the  least  response  to  wind. 

Closed-loop  Bank  Command  Frequency  Response 
The  frequency  responses  Y/Yrn  given  in  Figure  4-3  would  be  an  accu¬ 
rate  indication  of  system  radar  noise  response  if  the  system  were  linear. 
However,  the  effects  of  the  nonlinearities  in  the  aircraft  do  not  appear 
on  the  frequency-response  plots.  Perhaps  a  better  indication  of  the  ef¬ 
fects  of  radar  noise  on  the  system  response  is  given  in  the  frequency 
response  4>/Yrn,  where  4>  is  the  bank  command  and  Yrn  is  the  radar  noise 
signal.  The  high  frequency  signals  in  $  are  filtered  out  before  reaching 
Y,  the  aircraft  position.  However,  these  signals  can  cause  both  the 
rudder  and  the  ailerons  servos  to  reach  the  mechanical  limits  of  travel. 
This  limiting  is  a  nonlinear  effect. 

Shown  in  Figure  4-5  are  the  frequency  responses  for  five  of  the  six 
systems.  The  ydot  system  is  not  included  here.  Note  that  the  frequency 
responses  of  the  observer  systems  are  considerably  less  than  that  of  the 
a-B  system,  but  much  greater  than  that  of  the  half-gain  system. 

Conclusions:  If  the  mechanical  limiting  of  the  travel  of  the  rudder 
and  the  ailerons  is  a  problem,  then  the  observer  systems  should  perform 
better  than  the  a-B  system  with  respect  to  radar  noise  response. 

Monte  Carlo  Simulation  Results 


The  Monte  Carlo  simulation  discussed  in  [1]  was  used  to  determine 
actual  time  responses  to  radar  noise  and  to  wind.  The  wind  model  used 
was  simply  white  noise;  thus  the  use  of  a  more  realistic  wind  model  may 


change  the  results  somewhat.  The  radar  noise  model  used  was  that  given 
in  [3].  Twenty  simulations  were  run  for  each  Monte  Carlo  test,  and  the 
results  given  are  the  statistical  averages  of  these  twenty  runs. 

Given  in  Table  4-2  is  the  root-mean-square  values  of  aircraft  posi¬ 
tion  in  feet,  with  radar  noise  and  no  wind.  The  aircraft  flight  began 
at  a  range  of  13,223  feet  (60  seconds  flight  time  to  touchdown)  with  the 
aircraft  on  the  extended  centerline  of  the  runway;  i.e.,  y ( 0)  *  0.  Thus 
all  motion  of  the  aircraft  is  due  to  the  radar  noise,  and  hence  is  unde¬ 
sirable.  The  a-6  system  exhibits  the  largest  motion,  the  ydot  system 
the  smallest.  The  observer  systems’  motions  are  greater  than  that  of  the 
half-gain  system,  but  considerably  less  than  that  of  the  a-6  system. 

Table  4-3  gives  the  root-mean-square  values  of  the  bank  command  for 
the  same  simulations.  Note  that  the  data  follows  the  same  trends  as  that 
in  Table  4-2. 

Table  4-4  gives  additional  data  from  the  same  Monte  Carlo  simulation. 

This  table  gives  the  percentage  of  time  that  the  rudder  and  the  ailerons 

are  limited.  Note  that  limiting  does  not  occur  for  the  ydot  system,  and 

very  little  for  the  (z-0.85)  observer  system.  The  limiting  is  practi- 
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cally  the  same  for  the  half-gain  system  and  the  (z-0.8)  observer  system. 
Extensive  limiting  occurs  in  the  a-6  system. 

Table  4-5  gives  results  for  the  same  initial  conditions  on  the  air¬ 
craft,  except  for  this  simulation  the  radar  noise  signal  is  set  to  zero 
and  the  wind  input  is  white  noise  as  given  in  [1].  Note  that  the  a-6 
system  exhibits  the  smallest  motion  to  wind,  while  the  half-gain  system 
exhibits  the  largest.  The  observer  systems'  motion  falls  between  these 


TABLE  4-2 


Radar  Noise  Response-Position  Y 


system 

rms  aircraft 
position,  feet 

%  above 

half-qain  system 

cx-S 

2.52 

77 

half-gain 

1.42 

0 

ydot 

1.09 

- 

(z-0.85)3 

1.61 

13 

(z-0.8)3 

1.68 

18 

(z-0.75)3 

1.97 

39 

i 

1 
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TABLE  4-3 


Radar  Noise  Response-Bank  Command  $ 


sys  tem 

rms  bank 
command -degrees 

%  above 

ha If -gain  system 

a -B 

3.00 

104 

ha If -gain 

1.47 

0 

ydot 

0.17 

- 

(z-0.85)3 

1.60 

9% 

(z-0.8)3 

2.23 

52% 

(z-0.75)3 

2.70 

84% 
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TABLE  4-4 

Per  Cent  Time  Limited 

aileron 

rudder 

system 

limiting,  % 

1 imitinq 

a-B 

23 

42 

half-gain 

0 

16 

ydot 

0 

0 

(z-0.85)3 

0 

1 

(z-0.8)3 

2 

15 

(z-0.75)3 

9 

30 

I 

I 

i 

i 

5 

| 

"  * 

.1 

] 

4-14  J 
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TABLE  4-5 


Wind  Response-Position  Y 


system 

rms  aircraft 
position -feet 

%  above 
cx-B  system 

a-0 

1.90 

0 

ha If -gain 

3.35 

82 

ydot 

2.02 

6 

(2-0. 85)3 

2.36 

24 

(2-0.8)3 

2.16 

14 

(z-0.75)3 

2.05 

8 
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Conclusions:  The  (z-0.75)  observer  system  exhibits  excessive  mot  jn 
to  radar  noise,  and  the  (z-0.85)^  observer  system  exhibits  excessive  mo¬ 
tion  to  wind.  The  radar-noise  motion  of  the  (z-0.8)^  observer  system  ap¬ 
proaches  that  of  the  half-gain  system;  the  wind  motion  of  this  observer 

system  approaches  that  of  the  a-s  system.  In  addition,  the  initial-con- 
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dition  response  of  the  (z-0.8)  observer  system  approaches  that  of  the 
a-B  system. 


V.  MODEL  2  SYSTEM 


In  Chapter  4  an  extensive  comparison  was  made  between  six  different 
systems.  A  different  system  was  used  to  make  the  actual  closed-loop 
landings,  and  this  system  is  called  the  model  2  system  in  [3].  This  sys¬ 
tem  differs  significantly  from  the  a-B  system  and  the  half-gain  system 
given  in  Chapter  2,  and  will  now  be  described. 

The  block  diagram  of  the  a-B  filter-controller  of  the  a-8  system  is 
shown  in  Figure  5-1.  The  gains  ,  G^,  G^  and  G^  have  the  standard  termi¬ 
nology  given  in  Table  5-1.  The  differences  in  the  a-B  system  parameters 
and  the  model  2  system  parameters  are  given  in  Table  5-2.  Note  that  setting 
the  parameter  ^  to  unity  gives  that  a  filter  in  Figure  5-1  a  gain  of  unity; 
i.e.,  that  a  filter  is  removed.  Note  that  the  model  2  gains  in  the  deriva¬ 
tive  path  and  in  the  acceleration  path  have  been  reduced  from  those  of  the 
a-B  system.  Thus  we  expect  the  model  2  system  to  have  a  better  radar  noise 
response,  but  a  degraded  wind  response,  when  compared  to  the  a-e  system. 
However,  the  low-frequency  gain  of  the  model  2  system  is  unchanged  from  that 
of  the  a-B  system. 

In  the  following  developments,  the  model  2  system  will  be  compared  to  an 

observer  system.  However,  in  this  chapter,  the  observer  system  will  have  the 

same  controller  parameters  as  the  model  2  system;  i.e.,  the  observer  system 

Is  the  same  as  the  model  2  system,  except  that  the  a-B  filter  in  the  model 

2  system  is  replaced  with  an  observer.  Three  observers  will  once  again  be 
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considered,  with  characteristic  equations  of  either  (z-0.85)  ,  or  (z-0.8)  , 
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TABLE  5-1 


Gains  of  Figure  5-1 
_  K4L*KCL 

V  TIL 


G2  =  K1L*KCL 


g3  =  k2l*trl 


g4  =  k3l*tal 


TABLE  5-2 

Controller  Parameters 
(T  =  0.1  seconds) 


parameter 

a-S  system 
of  Chapter  2 

a-B  system 
model  2 

kcl 

0.1 

0.1 

TIL 

30.0 

30.0 

TRL 

7.5 

6.0 

TAL 

7.5 

2.5 

al 

0.3174 

1.0 

°2 

0.234 

0.234 

a3 

0.234 

0.234 

a4 

0.1211 

0.1211 

0.G1 

0.51 

B 

0.1746 

0.1746 

or  (z-0.75)' 


Initial  Condition  Response 

First  the  initial  condition  responses  will  be  compared.  For  the 
comparison,  both  the  radar  noise  input  and  the  wind  input  have  been  set 
to  zero. 

The  initial  condition  responses  are  given  in  Figure  5-2.  The  initial 
aircraft  displacement  y(0)  was  set  to  twenty  feet  at  a  range  of  13,223 
feet,  which  required  a  flight  time  to  touchdown  of  sixty  seconds.  The 
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response  of  the  (z-0.75)  observer  is  not  plotted,  since  it  is  very  close 
to  that  of  the  model  2  system. 

Conclusions:  The  (z-0.75)3  and  the  (z-0.8)3  observer  systems  degrade 

the  initial-condition  response  only  slightly.  Noticeable  degradation 
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occurs  in  the  (z-0.85)J  observer  syste.  .  Recall  that  the  (z-0.35)  obser¬ 
ver  is  the  slowest  one. 


Stability  Margins 

The  gain  margins  and  the  phase  margins  of  the  four  systems  are  given 
in  Table  5-3.  The  description  of  the  system  opened  at  different  points 
is  given  in  Chapter  4. 

Conclusions:  The  phase  margins  of  the  observer  systems  opened  at  <$>(t) 
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are  somewhat  small,  but  those  of  the  (z-0.8)  and  the  (z-0.75)  observer 
are  close  to  that  of  the  model  2  system  opened  at  $(k). 

Closed-loop  Noise  Frequency  Response 


Shown  in  Figure  5-3  are  the  closed-loop  frequency  responses  Y/Yrn> 
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TABLE  5-3 


system 
model  2 
(z-0.85)3 
(z-0.8)3 
(z-0.75)3 


(z-0.85)3 

(z-0.8)3 

(z-0.75)3 


Stability  Margins 
Opened  at  <*> ( k ) 

phase  margin 
37° 

41° 

41° 

41° 


gain  margin 
20  dB 
17  dB 

17  dB 

18  dB 


Opened  at  <j>(t) 


32° 

12 

dB 

35° 

15 

dB 

37° 

18 

dB 
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where  Y  is  the  aircraft  position  and  Y^  is  the  radar  noise  signal  (see 

rn 

3 

Figure  2-2).  The  frequency  response  of  the  (z-0.75)  observer  system 
is  not  plotted,  since  it  is  quite  close  to  that  of  the  model  2  system. 

Conclusions:  The  observer  systems'  noise  responses  degrade  from  that 
of  the  model  2  system,  with  the  slowest  observer  having  the  narrowest 
bandwidth. 


Closed-loop  Hind  Frequency  Response 

Given  in  Figure  5-4  is  the  closed-loop  frequency  response  Y/W,  where 
Y  is  the  aircraft  position  and  W  is  the  wind  input. 

Conclusions:  Since  low  wind  response  depends  on  a  wide  open-loop 
bandwidth,  we  expect  the  (z-0.85)^  observer  (slowest)  to  have  the  greatest 
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wind  response.  This  is  seen  to  be  true.  The  (z-0.75)  wind  response  is 
very  close  to  that  of  the  model  2  system,  and  is  not  plotted. 

Closed-loop  Bank  Command  Frequency  Response 

Given  in  Figure  5-5  are  the  frequency  responses  4>/Yrn  ,  where  $  is 
the  bank  command  and  Yrn  is  the  radar  noise  signal.  The  importance  of 
this  frequency  response  is  discussed  in  Chapter  4. 

Conclusions:  The  observer  systems  show  significant  improvement  in 
the  system  gain  from  radar  noise  to  bank  command.  Thus  the  ailerons  and 
rudder  should  exhibit  less  motion  in  the  observer  systems. 

Monte  Carlo  Simulation  Results 

The  Monte  Carlo  simulation  is  discussed  in  [1]  and  in  Chapter  4. 

Given  in  Table  5-4  are  the  root-mean-square  values  of  several  signals  obtained 
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from  Monte  Carlo  simulations  with  radar-noise  input.  The  column  labeled 
y  is  the  rms  value  of  position,  $  is  the  rms  value  of  bank  command,  s  is 
the  rms  aileron  angle,  is  the  rms  aileron  velocity,  5r  is  the  rms 
rudder  angle,  and  the  rms  rudder  velocity.  Also  given  are  the  per  cent 
of  time  that  the  aileron  servo  and  the  rudder  servo  are  limited.  Note 
that  the  observer  systems  show  significantly  less  response  to  radar  noise 
than  does  the  model  2  system. 

In  terms  of  radar-noise  response,  replacing  the  a- 8  filter  with  an 
observer  significantly  improves  response.  For  example,  the  (z-0.8)  ob¬ 
server  has  less  than  one-half  the  aileron  displacement  and  velocity  as  does 
the  model  2  system.  In  addition,  the  aileron  servo  limiting  for  this 
observer  system  is  almost  zero,  compared  to  13.3%  of  the  time  for  the  model 
2  system.  Figure  5-6  shows  a  typical  aileron  response  to  radar  noise  for 
the  model  2  system  and  for  the  (z-0.8)3  observer  system.  The  flight  time 
runs  from  30  to  46  seconds,  or  for  a  range  from  approximately  6600  feet  to 
3000  feet.  Note  that  the  aileron  motion  for  the  observer  system  is  signifi¬ 
cantly  less  than  that  for  the  model  2  system. 

Table  5-5  gives  results  from  Monte  Carlo  simulations  with  the  radar 
noise  removed  and  the  wind  input  as  given  in  [1],  Note  that  the  model  2 
system  exhibits  the  least  response  to  wind,  but  that  the  observer  systems' 
responses  are  not  significantly  greater.  Also  given  is  a  comparison  to 
the  wind  response  of  the  a-B  system  of  Chapter  4. 

Conclusions:  The  observer  systems  exhibit  a  significantly  improved 
radar-noise  response  when  compared  to  the  model  2  systems,  but  the  wind 
response  is  somewhat  degraded.  Thus,  if  an  observer  is  employed  in  the 
system,  probably  the  controller  gains  can  be  increased  to  reduce  wind  response 
while  maintaining  a  reasonable  radar-noise  response. 
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VI.  CONCLUSIONS 


The  F4J  aircraft  lateral  control  system  in  a  MATCALS  configuration  is 
investigated  in  this  report.  Several  different  controllers  are  utilized 
in  these  studies,  to  determine  which  yield  the  best  radar-noise  response 
and  which  yield  the  best  wind  response. 

The  proposed  MATCALS  systems  contain  an  a-8  filter  in  the  controller. 
Alternative  systems  are  constructed  by  replacing  the  a-8  filter  with  an 
observer. 

In  general  the  observer  control  systems  exhibit  significantly  less 
radar-noise  response  than  do  the  a-8  systems,  but  exhibit  somewhat  more 
wind  response.  These  studies  indicate  that  the  observer  controllers  im¬ 
prove  the  MATCALS  system's  operation  when  compared  to  the  a-8  controllers, 
and  that  the  observer  systems  should  be  considered  further. 
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THE  DESIGN  OF  A  TRI-STATE  ADAPTIVE  TRACKING 
FILTER  FOR  THE  MATCALS  SYSTEM 

ABSTRACT 

A  tri-state  adaptive  tracking  filter  is  designed  for  use  in  the 
F4J  aircraft  lateral  control  system  in  an  automatic  landing  configuration. 
The  system  presently  uses  an  alpha-beta  tracking  filter  to  estimate  the 
aircraft's  lateral  position  and  velocity.  The  tri-state  adaptive  filter 
is  designed  to  replace  the  alpha-beta  filter.  Results  from  the  F4J  lat¬ 
eral  control  system  simulation  •’’ndicates  that  the  tri -state  adaptive  fil¬ 
ter  improves  the  system's  response. 
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I.  INTRODUCTION 


This  study  investigates  the  possible  use  of  a  tri-state  adaptive 
tracking  filter  in  an  automatic  landing  system  for  aircraft.  The  auto¬ 
matic  landing  system  under  consideration  is  the  Marine  Air  Traffic  Con¬ 
trol  and  Landing  System  (MATCALS).  Presently  the  MATCALS  control  system 
uses  an  alpha-beta  filter  to  estimate  the  aircraft's  position  and  velo¬ 
city.  The  performance  of  this  system  is  greatly  degraded  by  the  amount 
of  noise  present  in  the  system.  The  tri-state  adaptive  tracking  filter 
is  presented  as  a  possible  replacement  for  the  alpha-beta  filter.  The 
performance  of  the  tri-state  adaptive  filter  was  evaluated  using  the 
simulation  of  the  F4J  aircraft  lateral  control  system. 

A  general  discussion  of  the  alpha  filter,  the  alpha-beta  filter, 
and  the  alpha-beta-gamma  filter  is  given  in  Chapter  II.  This  discussion 
includes  derivations,  stability  analysis,  and  noise  response  evaluation 
for  each  filter.  The  tri-state  adaptive  filter,  which  is  composed  of 
the  three  filters  discussed  in  Chapter  II,  is  presented  in  Chapter  III. 
The  selection  of  the  appropriate  filter  output  is  governed  by  the  vari¬ 
ance  of  the  filters'  smoothed  position  estimates. 

A  general  description  of  the  MATCALS  F40  Lateral  Control  System 
Simulation  is  given  in  Chapter  IV.  Chapter  V  presents  the  results  of 
the  F4J  lateral  control  system  simulation  using  the  tri-state  adaptive 
filter  in  place  of  the  alpha-beta  filter. 
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II.  THE  COMPONENT  FILTERS 


Introduction 

When  using  radar  track-while-scan  systems  to  track  a  moving  air¬ 
craft,  radar  returns  in  the  form  of  raw  digitized  position  measurements 
must  be  processed  to  provide 

1)  a  smoothed  estimate  of  present  position; 

2)  a  smoothed  estimate  of  present  velocity; 

3)  a  one  step  ahead  prediction  of  position  for  track 
correlation  or  bin  selection. 

In  the  design  of  these  tracking  systems,  two  conflicting  require¬ 
ments  must  be  met.  The  first  requirement  is  that  the  system  must  have 
good  noise  smoothing  properties.  A  system  of  this  type  is  typically 
characterized  by  a  sluggish  system  response,  long  time  constant  and  nar¬ 
row  bandwidth.  The  second  requirement  is  that  the  system  must  have  fast 
maneuver  following  capabilities.  This  type  of  system  is  usually  charac¬ 
terized  by  a  fast  system  response,  short  time  constant,  and  wide  band¬ 
width.  The  first  requirement  is  essential  because  of  the  inherent  noise 
present  in  the  raw  unprocessed  radar  position  measurements.  On  the 
other  hand,  in  tracking  airborne  targets,  the  dynamics  of  the  tracking 
system  must  be  capable  of  tracking  an  object  that  is  not  stationary  in 
space.  As  a  result,  the  second  requirement  is  necessary.  Unfortunately, 
as  one  makes  the  tracking  system  more  insensitive  to  noise,  the  maneuver¬ 
following  capabilities  suffer.  Thus,  some  compromise  is  always  required. 
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However,  the  smoothing  equations  should  be  structured  so  as  to  give  the 
best  compromise.  That  is,  noise  smoothing  should  be  maximum  for  a  given 
maneuver  following  capability  and  vice-versa. 

The  smoothing  of  the  position  and  velocity  estimates  is  usually  ac¬ 
complished  with  a  digital  filter.  These  filters  can  range  in  complexity 
from  a  simple  two  point  extrapolator  to  the  comparatively  complex  Kalman 
filter  [1].  This  section  examines  the  alpha  filter  [2],  alpha-beta 
filter  [3],  and  the  al pha-beta-gamma  filter  [4].  Although  the  Kalman 
filter  is  more  sophisticated  and  more  accurate  than  the  filters  being  ex¬ 
amined  in  this  section,  it  is,  computationally,  the  most  costly  to  imple¬ 
ment.  The  principle  advantage  of  the  alpha,  alpha-beta,  and  alpha-beta- 
gamma  filters  versus  the  Kalman  filter  is  that  the  computations  require 
no  matrix  inversions  as  does  the  Kalman  filter,  which  results  in  faster 
filter  outputs  and  less  data  storage. 

Throughout  this  discussion,  only  a  single  dimensional  filter  will 
be  considered.  The  variable  x  may  denote  the  aircraft's  position  in 
azimuth,  elevation,  or  range.  By  assuming  smoothing  in  any  given  coor-  : 
dinate  can  be  handled  independently,  extension  to  multiple  coordinate 
systems  is  straightforward. 

The  Alpha  Tracking  Filter 
Derivation  of  the  Alpha  Filter 

The  first  tracking  filter  to  be  examined  is  the  alpha  tracker. 

The  alpha  tracker  is  founded  upon  the  fundamental  theorem  of  exponential 
smoothing.  The  fundamental  theorem  states  that,  given  a  time  series 
{ X ( t ) }  with  observations  at  equally  spaced  intervals,  the  first  N+l  de¬ 
grees  of  exponential  smoothing  can  be  combined,  using  the  binomial 


*-  -  -—*■■:** V 
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coefficients,  to  give  an  estimate  of  the  values  of  the  coefficients  of  an 
Nth  degree  polynomial  model  of  the  observations  to  date,  evaluated  at  the 
time  of  the  most  recent  observation  [2].  The  alpha  tracker  takes  the 
form  of  the  simplest  first  degree  of  exponential  smoothing  and  can  be  ex¬ 
pressed  as  ' 

Xs(k)  =  a  xm(k)  +  (1-a)  x$(k-l )  .  (2-1)  1 

where, 

xs (k )  =  smoothed  estimate  at  time  kT 

xm(k)  =  input  measurement  of  position  at  time  kT 

a  =  position  smoothing  constant.  0  <  a  <  1  . 

Suppose  that  the  observations  of  position  to  date  give  an  estimate 
of  xs(k-l)  and  a  new  observation  x^Ck)  is  received  at  time  kT.  If  the 
new  value  is  larger  than  the  old  estimate,  the  new  estimate  x  (k)  will  be  * 

i  * 

I 

larger,  and  conversely. 

The  alpha  tracker  is  derived  from  the  zero  order  "window"  tracker  j 

[5],  Assume  that  there  exists  a  series  of  regularly  spaced  samples  of 
position  which  must  be  smoothed.  ! 

xm(k)’  xm(k_1)’  xm(k_2)’  •••»  xm(1>’  xm{0)  (2_2)  J 

such  that 

xm(k)  =  x(k)  +  n(k) 

where 

t 

xm(k)  =  input  measurement  of  position  at  time  kT 
x(k)  =  actual  value  (constant)  of  position  at  time  kT 
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n(k)  =  measurement  noise 

At  this  point  in  the  discussion,  the  noise  is  assumed  to  be  Gaussian, 

2 

zero  mean,  with  variance  =  a  .  The  noise  is  also  assumed  to  be  uncor- 

n 

related  from  sample  to  sample  and  therefore  independent  from  sample  to 
sample. 

If  x(k)  equals  some  constant  X,  then  the  "window  tracker"  with  a 

smoothing  interval  of  L  samples  might  be  utilized  to  "average  out"  the 

effects  of  the  random  noise  variations  on  x  (k)  such  that  the  estimate 

m 

is  approximately  equal  to  the  actual  value  X.  This  may  be  accomplished 
via  equation  (2-3). 

!  L 

x  (L)  =  f  I  x  (k)  (2-3) 

s  L  k=l  m 

Now,  if  the  approach  of  minimizing  mean-square  error  between  xs(L) 
and  xm(k)  is  taken,  the  best  least  squares  fit  of  a  horizontal  line  to 
the  data  is  obtained. 


min  {  z  (x  (L)  -  x _ ( i ) ) 2 }  =  0  (2-4) 

4  -i  s  m 

*s<u 


(xs(L)  -  xjl))2}  = 


L 

Z 

i=l 


2(x$(L)  -  xm(i))  =  0 

(2-5) 


L  L 

2  Z  x  (L)  -  2  z  x  (i )  =  0  (2-6) 

1=1  s  i=l  m 


L  L 

Z  x  (L)  =  z 
1=1  S  1=1 


(1) 


(2-7) 
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L  x  (L)  =  l  x  (1) 
s  i=l  m 


(2-8) 


xs(L)  *  f  xm(1> 


(2-9) 


Therefore,  the  "window"  tracker  provides  a  smoothed  position  estimate 
based  upon  the  minimization  of  the  mean-square  error  between  x$(L)  and 
xm(k). 

A  one  step  ahead  prediction  of  position  may  be  obtained  from  the 
"window"  tracker  by  merely  shifting  the  L-length  window  such  that  it 


spans  the  measured  data  from  xm(2)  to  xm(L+l).  Equation  (2-3)  now  takes 
the  form  of 


L+l 


VL*,)  =  t  ,f2x»(1) 


(2-10) 


L+l 


xpa*D  -  (i  ^  xra(D)  -  f  xm(i) 


(2-11) 


xp(L*l)  =  (i  ^  xm(l))  -  f  xmO )  +  f  xro(L-0) 


(2-12) 


xp(Ul)-xs(L)*J-(xm(L+l)-xB(l)) 


(2-13) 


It  is  safe  to  assume  that  the  smoothed  estimate  x$(L)  is  a  more  accurate 
value  of  true  position  than  the  noise  corrupted  measurement  xmO)- 
Therefore,  (2-13)  becomes 

x  (L+l)  =  x$(L)  +  1  (xm(L+l)  -  x$(D) 


(2-14) 
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rearranging  terms 

xp(L+1)  =  f  xm(L+1)  +  0  -  f}  XS(L)  (2_15) 

By  defining  <*  =  •£■  then  a  familiar  form  of  the  alpha  tracking  filter  is 
obtained  [6]. 

xp(k+l)  =  a  xm(k+1)  +  (1-a)  x$(k)  (2-16) 

where  x  (k+1)  is  the  one  step  ahead  prediction  of  position.  By  examining 
(2-15),  it  is  evident  that  as  the  length  L  of  the  smoothing  window  is 
extended,  a  decreases,  resulting  in  increased  smoothing  since  the  new 
measurement  data  is  not  weighted  as  heavily  as  the  previous  estimates  of 
position.  If  it  were  known  that  the  initial  estimates  of  positions  were 
correct  then  there  would  be  no  point  in  using  the  measurement  data  as  a 
basis  for  prediction  and  the  smoothing  constant  a  would  approach  the 
value  of  zero.  If,  on  the  other  hand  L  is  decreased,  a  increases  in 
value  (approaching  unity)  such  that  the  new  measurement  of  position  is 
weighted  more  heavily,  rapidly  discounting  the  effects  of  the  previous 
estimates. 

When  (2-16)  is  substituted  into  itself  for  successively  earlier 
samples,  a  series  emerges  of  the  form 

k+1 

x  (k+1 )  =  a  z  (1 -a)n  x  (k+1 -n)  (2-17) 

p  n=0  m 

From  (2-17)  it  is  evident  that  the  alpha  tracker  assigns  geometrically 
decreasing  weights  to  the  least  recent  data,  hence  the  term  exponential 
smoothing. 
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Alpha-Filter  Noise  Response 

The  alpha  tracker  can  be  depicted  as  in  Figure  2.1,  where  the  im¬ 
pulse  response  of  the  digital  filter  is  defined  by 

h(k)  £  a  (l-a)k  .  (2-18) 

The  output  of  the  filter  is  a  convolution  of  the  input  and  the  filter's 
impulse  response 

xs(k)  =  h(k)  *  xm(k)  .  (2-19) 

Some  idea  of  the  noise  smoothing  performance  of  the  alpha  tracker  may  be 
obtained  by  examining  the  mean  and  variance  of  xs(k).  The  convolution  of 
(2-19)  may  be  written  as 


x  (k)  =  a  z  (l-a)n  [X  +  n(k-n)].  (2-20) 

s  n=0 


where 

a(l-a)n  =  impulse  response  of  alpha  tracker 

X  =  actual  position  of  target 

n(k)  =  additive  zero  mean  Gaussian  noise  with 

variance  =  a  2 
n 

Taking  the  expected  value  of  (2-20)  yields  the  mean  of  the  output 
of  the  alpha  tracker. 


E{x  (k) }  =  E{a  I  (l-a)n  [X  +  n(k-n)]}  (2-21) 

s  n=0 


E{x  (k) }  =  a  Z  (l-a)n  E{X  +  n(k-n)} 
s  n=0 


(2-22) 
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E{x  (k) }  =  a  E  ( 1 -a ) '  E { X }  +  E{n(k-n) } 
s  n=0 


(2-23) 


Since  X  is  a  constant,  E { X }  =  X.  n(k)  is  zero  mean,  therefore  E{n(k-n) > 
=  0  and  (2-23)  becomes 


E(x  (k) }  =  a  I  (1 -a)n  X 
s  n=0 


(2-24) 


Using  the  identity,  for  j a !  <  1, 


n  1 


1  -a 


(2-25) 


E{x$(k)}  = 


X  =  X 


(2-26) 


Thus  the  alpha  tracker  is  capable  of  tracking  a  constant  level  of  posi¬ 
tion,  corrupted  by  zero  mean  Gaussian  noise,  with  no  bias  or  offset. 

The  variance  of  the  output  estimate  x£(k)  is  given  by 


var{x  (k) }  =  E{[a  i  (l-a)n  [X+n(k-n)]  -  E{xc(k)}f}  .  (2-27) 

s  n=0  s 

From  (2-26)  E{xs (k) }  =  X 
so  that 


OO  A 

varix  (k ) }  =  E{[a  z  (l-a)n  (X+n(k-n)  -  X)]2} 
5  n=0 


(2-28) 


var{x  (k)}  =  Ef[a  z  (l-a)n  n(k-n)]2} 
s  n=0 


(2-29) 


var{x  (k) }  =  a2  E{[  Z  (l-a)n  n(k-n)][  z  (1-a)^  n(k-j)]}  (2-30) 
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var{x  (k) }  =  2  z  l  (l-a)n  (l-a)j  E{n (k-n )  n(k-j)}  (2-31) 

n=0  j=0 

The  noise  n(k)  is  uncorrelated  from  sample  to  sample  such  that 

E{n(k-n)  n(k-j)}  =  0  if  i  t  j 
E{ n ( k-n )  n(k-j)}  =  o 2  if  i  =  j 

Therefore  (2-31)  is  reduced  to  only  one  summation. 

var{x  (k) }  =  a2  l  (l-a)2n  a  2  •  (2-32) 

s  n=0  n 

“  k  1 

Using  the  identity  l  a  =  —  again, 
n=0  l_a 

2 


var{x  (k)}  =  - 2 - 

(2-33) 

S  l-(l-a)2 

var{x  (k)}  =  — -  c  2 

S  c-a  n 

(2-34) 

It  is  evident  that  as  a  decreases  the  noise  response  of  the  system,  and 
thus  the  variance  in  the  output  estimate,  decreases. 

Analysis  of  the  Alpha  Filter  in 
the  Z-aomain 

By  taking  the  z-transform  [7]  of  equation  (2-1)  the  alpha-tracker 
may  be  examined  in  the  z-domain, 

x$(z)  =  (l-a)z"1  x$(z)  +  a  xm(z) 


(2-35) 
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Rearranging  terms  yields  the  z-transform  of  the  alpha-tracker  [6]. 


X  (2) 

H(z)  =  irjT) 

m 


_ a _ 

1  - ( 1 -a )z-1 


The  alpha-tracker  is  seen  to  have  a  pole  located  at 
z  =  1-a 


(2-36) 


(2-37) 


Since  a  tracking  system  with  stable  dynamics  is  required,  a  must  be  se¬ 
lected  such  that  the  system  pole  lies  within  the  unit  circle.  The  locus 
of  possible  pole  locations  within  the  unit  circle  is  shown  in  Figure  2.2. 
By  examining  Figure  2.2  it  is  evident  that  a  must  fall  within  the  range 
of  0  <  a  <  1  for  stable  operation.  Selecting  a  too  close  to  1  causes  the 
weighting  sequence  of  the  alpha-tracker  to  rapidly  approach  zero.  Such  a 
selection  yields  a  system  with  fast  response  to  transients.  Selecting  a 
close  to  0  results  in  a  sluggish  system  response  (more  smoothing)  and  a 
decreased  response  to  noise.  The  zero  of  the  alpha-tracker  is  located  at 

z  =  0  (2-38) 

and  thus  weights  the  entire  frequency  response  of  the  system  equally. 

The  alpha-tracker,  though  simple  to  implement,  is  limited  in  its 
ability  to  track  a  moving  target.  The  limited  tracking  ability  of  the 
alpha-tracker  becomes  evident  when  xs(k)  is  examined  as  a  function  of 
xm(k)  for  different  target  trajectories . 

The  case  of  constant  position,  zero  velocity  will  first  be  exam¬ 
ined.  This  case  might  be  analogous  to  the  azimuth  or  elevation  tracking 
of  an  aircraft  moving  radially  with  respect  to  the  radar. 


xm(k>  =  X 


(2-39) 


Figure  2.2.  Locus  of  Possible  Pole  Locations  for  Stable 
Operation  of  Alpha  Tracking  Filter. 
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Applying  the  z  transform, 

_  Xz 

Xn/Z^  "  z-1 

(2-40) 

"W  -rffer  • 

(2-41) 

Therefore 

/ ~ \  XZ  aZ  _  aXz 

XSU;  "  Z-1  '  z-(l-a)  (z-l)(z-(l-ct))  • 

(2-42) 

By  utilizing  the  partial  function  expansion  method, 

x  f2)  ,  Xz  _  (l-g)Xz 
xsU)  z-1  z-(l-o)  ‘ 

(2-43) 

Finally, 

Xs(k)  =  X  -  (l-a)(l-a)kX 

(2-44) 

Xs(k)  =  X  -  (l-a)k+1  X  , 

(2-45) 

where 

X  =  actual  signal 

k+1 

X(l-a)  =  bias  or  amount  of  lag  . 

As  k  increases,  the  bias  rapidly  approaches  zero  and  the  output  equals 

the  constant  position  input  X. 

The  alpha-tracker  does  not  perform  as  well  when  a  constant  velocity 

input  is  applied  to  the  system.  This  can  be  verified  by  making  Xfn(k)  a 

sampled  ramp  function. 

x  (k)  =  Vk 
m 


(2-46) 
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XJZ) 


Vz 


(z-l)‘ 


Therefore, 


x  (z)  =  - - 

S  (z-1  )2 (z- ( 1 -a )  ) 


After  applying  the  partial  fraction  expansion, 


x  (k)  =  Vk 


V(l-g) 


[l-(l-aj] 


(2-47) 


(2-48) 


(2-49) 


where, 

Vk  =  actual  signal 

[Ml-ct)k]  «  bias  . 

Unlike  before,  as  k  increases,  the  bias  does  not  approach  zero  such  that 

x  (k)  =  Vk  -  VUral  .  (2-50) 

s  a 

From  (2-49)  the  inability  of  the  alpha-tracker  to  track,  with  zero  bias, 
any  signal  other  than  a  constant  amplitude  signal  with  random  noise  vari¬ 
ations  is  established. 


The  Alpha-Beta  Tracking  Filter 
Derivation  of  the  Alpha-Beta  Filter 

The  inability  of  the  alpha-tracker  to  track  without  bias  a  target 
whose  position  is  changing  with  respect  to  time  suggests  that  a  more 
generalized  model  for  target  trajectory  should  be  investigated.  The 
new  model  should  account  for  the  target's  change  in  position  via  a  con¬ 
stant  velocity  term  such  that  the  new  model  for  target  trajectory  is 
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x(k)  =  x(k-l )  +  vT  .  (2-51 ) 

where 

x(k)  =  actual  position  at  time  kT 

x(k-1)  =  actual  position  at  time  (k-l)T 

v  =  velocity  (constant) 

T  =  time  between  samples 

To  track  a  target  following  a  trajectory  as  prescribed  in  (2-51), 
both  smoothed  estimates  of  position  and  velocity  are  required.  From  [3] 
one  method  to  obtain  smoothed  estimates  of  position  and  velocity  is  to 
define  a  tracking  system  such  that 

xs(k)  =  gx(0)xm(k)  +  gx(l)xm(k-l)  +  ...  (2-52) 

and 

x$(k)  =  gx(0)xm(k)  +  gx(l)xm(k-l)  +  ...  (2-53) 

where 

x  (k)  =  noise  corrupted  position  measurement  at  time  nT 
(xm(k)  =  x(k)  +  n) 

xs(k)  =  smoothed  estimate  of  position  at  time  kT 
xs ( k )  =  smoothed  estimate  of  velocity  at  time  kT 
gx(k.)  =  the  input  position  to  output  position  weighting 
sequence  or  unit-impulse  response 
gx(k)  =  the  input  position  to  output  velocity  weighting 
sequence  or  unit-impulse  response 
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Smoothed  estimates  of  position  and  velocity  may  also  be  obtained  by 
defining  a  set  of  recursive  equations  of  the  form 

xs(k)  =  xs(k-1)  +  y2  x$ ( k-2)  +  ...  +  yn  xs ( k-N ) 

*  60  xm(k>  +  xn,(k-')  +  •••  +  5m  xn/k-",)  <2'54) 

x$(k)  =  n-j  xs(k-l)  +  n2  x$(k-2)  +  ...  +  nN  x$(k-N) 

♦  X0  xm(k)  *  X,  xm(k-1 )♦...+  Xm  xm(k-m)  .  (2-55) 

Equations  (2-54)  and  (2-55)  are  termed  Nth  order  (since  Yk  and  =  0  for 
k  >_  N+l )  and  at  least  N  storage  locations  in  a  computer  are  required  to 
compute  these  estimates.  Equations  (2-54)  and  (2-55)  may  be  combined  to 
form  a  set  of  second  order  difference  equations  known  as  the  alpha-beta 
tracking  filter  [8]. 

xs(k)  =  Xp(k)  +  a(xm(k)  -  Xp(k) )  (2-56) 

x$(k)  =  x$(k-l )  +  e/T  (xm(k)  -  Xp(k))  (2-57) 

Xp(k+1)  =  x$(k)  +  T  xg(k)  (2-58) 

xp(k+l )  =  x$(k)  (2-59) 

where 

xm(k)  =  noise  corrupted  position  measurement  at  time  kT 
x$(k)  =  smoothed  estimate  of  position  at  time  kT 

x$(k)  =  smoothed  estimate  of  velocity  at  time  kT 

x  (k+1)  -  one  step  ahead  prediction  of  position  at  time 


(k+1  )T 
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Xp(k+1)  =  one  step  ahead  prediction  of  velocity  at  time 
(k+l)T 

a  =  position  smoothing  constant 

6  =  velocity  smoothing  constant 

T  =  time  between  samples 

Successive  corrections  are  made  to  the  smoothed  values  of  position 
and  velocity  proportional  to  the  differences  between  the  measured  posi¬ 
tion  *m{k)  and  the  previous  predicted  estimate  of  position  xp(k).  As 
time  progresses  from  the  interval  (k-l)T  <  t  ^  kT  to  the  next  interval, 
kT  <  t  <_  (k+l)T,  x$(k)  is  increased  by  the  amount  a[xm(k)  -  xp(k)],  and 
the  velocity  xg(k)  is  increased  by  the  amount  3[xm(k)  -  xp(k)].  Examin¬ 
ing  (2-56)  through  (2-59)  it  is  apparent  that  the  value  of  a  and  3  deter¬ 
mine  the  degree  of  smoothing  between  measured  and  predicted  values  of 
position  to  yield  the  smoothed  estimate  of  position.  If  a  and  3  approach 
unity,  then  the  measured  values  of  position  are  weighted  more  heavily 
than  the  predicted  values  of  position.  Such  a  system  is  characterized 
by  wide  bandwidth  resulting  in  very  little  smoothing  of  the  incoming  po¬ 
sitional  data.  Conversely,  if  a  and  3  approach  zero  then  the  predicted 
values  of  position  are  weighted  more  heavily  than  the  measured  values 
effectively  narrowing  the  bandwidth  of  the  system  and  resulting  in  in¬ 
creased  data  smoothing.  The  predictor  equations  are  based  upon  the  as¬ 
sumption  that  the  target  aircraft  will  maintain  a  constant  velocity 
trajectory  throughout  the  tracking  procedure. 
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Analysis  of  the  Alpha-Beta  Filter 
in  the  Z-Domain 

•The  predictor  equation  (2-5d)  can  be  modeled  as  a  sampled 
data  feedback  system  in  which  the  error  defined  by  e  (t)  =  x  (t)  - 

A  III 

Xp(t),  is  converted  to  a  train  of  impulses, 

e*(t)  £  E  e  (kT)  6(t-kT)  .  (2-60) 

x  k  x 

This  train  is  applied  to  a  single  and  double  integrator  combination  whose 
impulse  response  is  the  sum  of  a  step  function  of  height  a  and  a  ramp  of 
slope  b/T.  The  block  diagram  is  shown  in  Figure  2.3. 

The  transfer  function  of  the  single  and  double  integrator  combina¬ 
tion  may  be  expressed  via  its  Laplace  transform  as 

Gn(s)=f  +  -^  ,  (2-61) 

The  transformation  of  the  continuous  error  ex(t)  into  the  impulse  train 
e*(t)  is  indicated  by  the  sampling  switch. 

In  the  z-domain,  the  transfer  function  of  the  predictor  can  be  ex¬ 
pressed  as 

G  (z)  =  I  C(a+Bk)z_k]  -  a  ,  (2-62) 

p  k=0 

which  suns  to  the  following  closed  form  [9], 

G  (z)  -  k*B>z~Y«~2  .  (2-63) 

p  o-z-'r 

By  manipulation  of  (2-63),  the  predictor  transfer  function  can  be  found 


to  be 


eedbac 

iha-Bet 
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H3(z) 

.  xp(2)  . 

(a+6 )z  -  a 

(2-64) 

V17 

z2  + 

(0+ a-2)z  +  (1-a) 

In  a  similar  manner,  the 

transfer  functions  for  smoothed 

position  and 

smoothed  velocity  may  be 

found 

to  be 

H2(z) 

Mz) 

0 

z(z-l ) 

(2-65) 

7z2 

+  (0+a-2)z  +  (1-a) 

H-|  (z) 

X$(z) 

z[aZ  +  0  -  a] 

(2-66) 

y77 

z2  + 

(0+a-2)z  +  (1-a) 

Note  that  the  three  transfer  functions  have  the  same  characteristic 
equation  and  the  poles  of  the  transfer  functions  can  be  found  solving 
for  the  roots  of, 

z 2  +  ( B+ot-2  )z  +  (1  -a)  =  0  .  (2-67) 

By  applying  the  quadratic  formula,  the  roots  of  the  characteristic  equa¬ 
tion  are  seen  to  occur  at 

Z  =  ~ =£  +  J  (b2  +  a2  +  2a0  -  4b)1/2  .  (2-68) 

Stability  of  the  Alpha-Beta  Filter 

Since  a  stable  tracking  system  is  desired,  the  values  for  a  and  0 
must  be  selected  such  that  the  system  poles  lie  within  the  unit  circle  of 
the  z-plane.  One  method  of  determining  the  values  of  a  and  0  for  stabil¬ 
ity  is  to  apply  the  Jury  stability  test  [10].  The  Jury  stability  test 
is  a  stability  criterion  for  sampled  data  systems  that  is  similar  to 
the  Routh-Hurwi tz  criterion  [11].  Applying  the  Jury  criterion  to  (2-67), 
the  Jury  array  obtained  is 
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0 

TT^T 


1 


z 

(ct+B-2) 


(2-69) 


The  necessary  and  sufficient  conditions  for  the  second  order  characteris¬ 
tic  equation  (2-67)  to  have  no  roots  on  or  outside  the  unit  circle  are  as 
fol lows : 


Condition  1 
Condition  2 
Condition  3 


F(1 )  >  0 
(-1 )2F(-1 )  >  0 


Condition  1 : 


z  +  (a+B-2)z  +  (1-ct) 


>  0 


z=l 


(1)^  +  (o+B-2)(l)  +  (1 -a)  >  0 

B  >  0 


(2-70) 


Condition  2:  (-1)2  [z2  +  (ct+3-2)z  +  (1  -a ) ] 

1  +  2-a+B  +  1 -a  >  0 
2a+B  <  4 


>  0 


z=-l 


I 

1 

l 


Condition  3:  ]a^|  < 

where. 


n-«i  <  1 


therefore 

a  >  0. 


An  additional  stable  condition  exists  when  6=0,  since  in  that  case  the 
zero  of  the  denominators  of  equations  (2-64)  through  (2-66)  at  z  =  1  is 
cancelled  by  one  of  the  zeroes  of  the  numerators.  The  resulting  require¬ 
ments  for  stability  are 

a  >  0 

6  >  0 

2a  +  0  >  4  .  (2-71) 


Alpha-Seta  Filter  Noise  Response 

Some  insight  into  the  noise  smoothing  performance  of  the  alpha-beta 
filter  may  be  obtained  by  examining  the  mean  and  variance  of  the  smoothed 
position  estimate  xs(k).  Utilizing  the  partial  fraction  expansion  method 
to  inverse  transform  equation  (2-66),  the  impulse  response  for  the 
smoothed  position  output  of  the  alpha-beta  filter  is. 


a  -aS+28 


2(0  +a+2a8-48) 


,  1  (B2.a2+a6-46)1/2:k 


*  If 


-a  -a6+2B _ 1  r-2-a-S 

2(s2+a2+2a0-4B)^22^  2 


1  (sV^e-os)122]1' 


(2-72) 


Assuming  a  constant  velocity  target,  the  input  to  the  filter  is  given  by: 


xm(k)  =  X  +  vkT  +  n(k) 


(2-73) 


where, 


X  =  actual  position  of  target  aircraft 
v  =  actual  velocity  of  target  aircraft 
n(k)  =  additive  zero  mean  Gaussian  noise  with  variance  =  o 
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T  =  sampling  interval  , 

The  smoothed  position  estimate  x$(k)  is  given  by  the  convolution  of  the 
impulse  response  h^(k)  and  the  input  xm(k). 


x  (k)  =  z  h, (n)  xm(k-n) 
n=0 


(2-74) 


From  (2-72)  the  impulse  response  may  be  expressed  as 


h-j  (k)  =  C1[B1]k  +  C2[B2]k 


(2-75) 


where, 


-a  -ctB+2g  _ 

2  (02+a2+2ae-4B)1/2 


-a  -aft+2s _ 

2  (f2+a2+2a8-4s )^  ^2 


2-ot-8  ,  1  - „2 .  2 


B]  =  ~  (B  +a  +2a6-4g)1 


B2  =  ~  -  \  (62+a2+2aB-4B)1/2 


Substituting  (2-75)  and  (2-73)  for  h^(k)  and  xm(k)  respectively,  the 
smoothed  position  estimate  is  given  by 


x  (k )  =  z  [C,(B1)n  +  C?(Bp)n][X  +  vkT  +  n(k-n)]  .  (2-76) 

s  n=0  c  c 

The  mean  of  x  (k)  may  be  found  by  taking  the  expected  value  of  (2-76) 


E{x  (k } }  =  E{  z  [C,(B, )n  +  C?(B?)n][X  +  vkT  +  n(k-n)]}  (2-77) 

s  n=0  ^ 
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E{x$(k)}  = 

E{xs(k) }  = 

Since  X  and  v  are  deterministic,  and  n(k)  is  zero  mean, 

E  { X  >  =  X 

E{vkT }  =  vkT 

E{n(k-n)}  =  0 

such  that  (2-79)  becomes 


£  [C-|  (B.j  )n  +  C?(B?)n]  E { [ X  +  vkT  +  n(k-n)]} 


n=0 


"2'  2 


(2-78) 


I  [C, (B1 )n  +  C2(B2)n][E{X)  +  E{vkT }  +  E{n(k-n)}] 
n=0 


ro  7o\ 


E{x$(k)} 
E  { x  s  ( k ) } 
Using  the  identity 
E(x$(k) } 


=  I  CC-j  (B-j  )n  +  C2(B2)n][X  +  vkT] 
n=0 


[X  +  vkT]  E  [C, (B, )n  +  C?(B?)n] 
n=0 


E  a 
n=0 


k  _  1 


1  -a  ’ 


=  [X  +  vkT] 


C]  +  C2 


.1_B1  1_B2. 


(2-80) 

(2-81) 

(2-82) 


After  extensive  algebraic  manipulation,  it  can  be  shown  that  the  mean  of 
the  smoothed  position  estimates  is 


E{x$(k)}  =  X  +  vkT 


(2-83). 


The  result  shown  in  (2-83)  implies  that,  given  a  constant  velocity  target 
trajectory,  the  -moothed  position  estimate  x  (k)  is  unbiased. 
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In  a  similar  manner,  the  variance  of  the  smoothed  position  estimate 
may  be  calculated. 


var{x  (k) }  *  E{[[  z  C,(B  )n  +  C?(B?)n][X  +  vkT  +  n(k-n) 
s  n=0 


-  E{xs(k)]]2} 


(2-84) 


From  (2-83) 


E { x  s ( k ) }  =  X  +  vkT  , 


so  that  (2-84)  becomes 


var{x$(k)} 


E{[  I  (C/B,)"  +  C?(B  )n)n(k-n)]2}  (2-85) 

n=0  c  ^ 


var{x$(k) } 


E{[Ci  z  (Bi  )nn(k-n)  +  C9  z  (Bjnn(k-n)] 
n=Q  c  n=0  c 


[C,  I  (B,)Jn(k-j)  +  C?  •  (B?)jn(k-j)]}  (2-86) 

'  j=0  1  j=0  c 


var{xs(k)} 


l  E  [C,(B,)n  +  C2(Bp)n][C1(B1)j 
n=0  j=0  ^  11 


+  Cp(Bp)J]  E{n(k-n)n(k-j)} 


(2-87) 


The  noise  is  uncorrelated  from  sample  to  sample,  therefore 

E{n(k-n)  n(k-j)}  =  0  if  i7j 
E{n(k-n)  n(k-j)}  =  op2  if  i=j 

Now  (2-87)  is  reduced  to  a  single  sum 


var{x  (k) }  =  X  [C. (8. )n  +  C?(B0)n]2  0„2 
s  n=Q  1  1  2c  n 


(2-88) 


I 


var{xs(k)}  =  £  [C2^)2"  +  2C]C2(B1  )n(B2)n  +  c|(B2)2n]  op2 


(2-89) 


Since  (B.)n(B0)n  =  (B,B0)n  and  z  ak  =  ~~  ,  (2-89)  becomes 
1  1  1  1  n=0  '"a 

'c  2  2C  C?  C?2  1  2 

var{Xs(k)) .  +  _  j  „n  _ 


(2-90) 


With  considerable  algebraic  manipulation,  it  can  be  shown  [3]  that 
variance  of  the  smoothed  position  estimate  is 


2:  (4^f J 


(2-91) 


If  the  variance  in  the  smoothed  position  is  expressed  as 


var{xs(k)}  =  Kx(0)an^ 


(2-92) 


then  Kx(0)  is  known  as  the  variance  reduction  ratio  [9],  If  K  (0)  is 
less  than  one  then  the  tracking  system  has  reduced  the  effect  of  noise 
on  the  output  signal.  If  ( 0)  is  greater  than  one,  the  system  has  ac¬ 
centuated  the  noise.  The  variance  reduction  ratio  is  superior  for  values 
of  a  and  6  approaching  zero. 


Performance  Measures  for  the 
Alpha-Beta  Filter 

Benedict  and  Border  [3]  introduced  two  performance  measures  in 
order  to  assess  properly  the  attributes  of  noise  reduction  and  transient 
performance  and  derive  an  optimal  relation  between  a  and  8. 

For  noise  smoothing,  the  performance  measures  are  the  aforemen¬ 
tioned  variance  reduction  ratio  in  position  output  defined  by: 
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«x(0)  =  steady-state  variance  in  position  output  , 

variance  in  raw  position  input  (2-93) 

and  the  variance  reduction  ratio  in  velocity  output  defined  by: 

K^(0)  =  steady-state  variance  in  velocity  output 

variance  in  raw  position  input  (2-94) 

For  transient  performance  the  performance  measures  are 


2  “  2 

D  =  T.  [(unit-increment  ramp )- (posi tion  ramp  response)] 

x  n=0 


(2-95) 


and 


2 

D-  =  Z  [(velocity  of  unit-increment  ramp) 
x  n=0 

2 

-  (velocity  ramp  response)]  .  (2-96) 

If  the  impulse  responses  of  both  the  smoothed  position  output  (previously 
derived  in  this  paper)  and  the  smoothed  velocity  output  are  obtained, 
then  the  performance  measures  for  the  alpha-beta  tracker  is  found  to 
be  [3], 


Kx(0> 

_  2a  +  s(2-3a) 
a[4-0-2a] 

(2-97) 

K-<0) 

_  1  2B2 

af4-2a-B] 

(2-98) 

II 

CM 

X 

o 

aB[4-B-2a] 

(2-99) 

II 

CM 

X 

o 

1  a2(2-a)  +  2s(l-a) 

(2-100) 

t2  a8t4-B-2a] 
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Based  upon  these  performance  measures,  Benedict  and  Bordner  [3] 

utilized  the  calculus  of  variations  to  optimally  select  the  values  of  a 

2  2 

and  g  so  as  to  minimize  respectively  Dx  and  for  a  given  Kx(0)  and 

K^(0),  and  vice  versa.  The  relationship  between  a  and  6  which  results 
in  the  optimal  tracker  is  given  by 


The  recursive  nature  of  the  alpha-beta  filter  equations  allows  them 
to  be  incorporated  easily  into  a  computer  algorithm  which  performs  the 
filtering  process.  A  flowchart  for  implementing  the  alpha-beta  filter  on 
a  computer  is  given  in  Figure  2.4. 

The  Alpha -Beta-Gamma  Tracking  Filter 

Derivation  of  the  Alpha-Beta- 
Gamma  Tracking  Filter 

The  alpha-beta  tracking  filter  previously  discussed  is  best  suited 
for  a  target  under  track  which  has  constant  velocity.  For  the  case  of 
an  aircraft  experiencing  a  great  deal  of  maneuvering,  changes  in  velocity 
between  sample  intervals  are  often  significant.  As  a  result,  the  ability 
of  the  alpha-beta  filter  to  track  an  accelerating  target  is  severely 
hindered.  The  degradation  in  performance  which  occurs  when  an  alpha-beta 
filter  tracks  an  accelerating  target,  suggests  that  a  tracking  filter 
which  includes  an  acceleration  estimate  might  be  appropriate  for  a  ma¬ 
neuvering  aircraft.  The  model  for  the  accelerating  target  trajectory  is 
therefore  chosen  as 

2 

x(k)  =  x  ( k- 1 )  +  vT  +  , 


I 


(2-102) 
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Figure  2.4.  Flowchart  of  the  Alpha-Beta  Tracking  Filter 
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where, 

x(k)  =  actual  target  position  at  time  kT 

x(k-l)  =  actual  target  position  at  time  (k-l)T 

v  =  velocity  of  target 

a  =  acceleration  of  target  (constant) 

T  =  time  between  samples 


In  order  to  accurately  track  a  target  following  a  trajectory  such 
as  prescribed  in  (2-102),  smoothed  estimates  of  position,  velocity  and 
acceleration  must  be  available.  Simpson  [4]  defines  a  system  of 
weighting  sequences  in  order  to  obtain  the  smoothed  estimates  of  posi¬ 
tion,  velocity,  and  acceleration.  These  are: 


where 


x  (k)  =  z  g  (n)  x  (k-n) 
s  n=0  x  m 


(2-103) 


x  (k)  =  z  g.(n)  x  (k-n) 
s  n=0 


(2-104) 


x" (k)  =  z  g--(n)  x  (k-n) 
s  n=0 


(2-105) 


*m<k> 


x$(k) 

*s(k> 

Vk> 

9X(") 


=  noise  corrupted  position  measurement  at  time  kT 
(xm(k)  =  x(k)  +n). 

=  smoothed  estimate  of  position  at  time  kT 
=  smoothe  r^timate  of  velocity  at  time  kT 
=  smoother  .timate  of  acceleration  at  time  kT 
=  the  input  position  to  output  position  weighting 
sequence,  or  unit-impulse  response  of  the  system 
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g^(n)  =  the  input  position  to  output  velocity  weighting 
sequence,  or  unit-impulse  response  of  the  system 
gx ' ( n )  =  the  input  position  to  output  acceleration  weighting 
sequence,  or  unit-impulse  response  of  the  system  . 

Alternately,  a  corresponding  sampled  data  system  can  be  character¬ 
ized  by  a  set  of  recursive  equations  such  as 

xs(k)  =  p1  xs(k-l)  +  p£  xs(k-2)  +  ...  pN  xs(k-N) 

+  <Sg  xm(k)  +  61  xm(k-l)  +  ...  +  6m  xm(k-m)  (2-106) 

x$(k)  =  n1  x s ( k- 1 )  +  n2  x$(k-2)  +  ...  +  nN  x$(k-N) 

+  \n  x  (k )  +  x  (k-1 )  +  ...  +  x  x  (k-m)  (2-107) 

0  m  1  m  mm 

x*s(k)  =  n1  x* s ( k-1 )  +  n2  x*s(k-2)  +  ...  +  nN  x's(k-N) 

+  e0  xm(k>  *  E1  xm(k-'>  +  £m  xn,(k-n,)  •  <2-,08> 

When  such  recursive  equations  terminate  (i.e.,  =  0  for 

k  N+l ) ,  the  equations  are  termed  Nth  order. 

Equations  (2-106)  through  (2-108)  may  be  combined  to  form  a  set  of 
third  order  difference  equations  which  define  the  alpha-beta-gamma  track¬ 


ing  filter  [4].  These  are  given  as 

x$(k)  =  xp(k)  +  a(xm(k)  -  xp(k))  (2-109) 

x$(k)  =  xp(k)  +  (s/T)  (xm(k)  -  xp(k))  (2-110) 

x's(k)  =  x*$(k-l )  +  (y/T2)  (xm(k)  -  xp(k) )  (2-111) 

xp(k)  -  x$(k-l)  +  T  x$(k-l)  (2-112) 
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x  (k)  =  xs(k-l)  +  T  x'(k-l)  ,  (2-113) 

where 

xm(k)  =  noise  corrupted  position  measurement  at  time  kT 
xs(k)  =  smoothed  estimate  of  position  at  time  kT 

x$(k)  =  smoothed  estimate  of  velocity  at  time  kT 

x’s(k)  =  smoothed  estimate  of  acceleration  at  time  kT 
Xp(k)  =  one-step  ahead  prediction  of  position  at  time  kT 

Xp(k)  =  one-step  ahead  prediction  of  velocity  at  time  kT 

a  =  position  smoothing  constant 

6  =  velocity  smoothing  constant 

y  =  acceleration  smoothing  constant 

T  =  time  between  samples 

The  alpha-beta-gamma  tracking  filter  is  a  third  order  filter  which 
can  track  a  constant  acceleration  target  with  a  trajectory  as  prescribed 
in  (2-102)  with  zero  steady-sta+e  error.  The  similarities  between  the 
alpha-beta  filter  and  the  alpha-beta-garrma  filter  are  obvious  when 
equations  (2-109)  through  (2-113)  are  compared  to  those  for  the  alpha- 
beta  filter.  While  a  and  6  perform  the  same  function  as  they  did  in  the 
alpha-beta  filter,  the  y  term  provides  the  much  needed  acceleration  esti¬ 
mate  essential  for  accurately  tracking  a  maneuvering  target.  This  more 
generalized  tracking  filter  not  only  maintains  track  throughout  maneuvers 
or  turbulent  conditions,  but  in  addition,  provides  smoothed  estimates  of 
position,  velocity,  and  acceleration  with  very  little  increase  in  com¬ 
putational  effort. 
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The  Alpha-Beta-Gamma  Filter 
in  the  Z- Domain 

The  properties  of  the  alpha-beta-gamma  tracking  filter  may  be  more 
effectively  studied  in  the  Z-domain.  Taking  the  Z-transform  of  equations 
(2-109)  through  (2-113)  and  using  the  right  shifting  property  results  in 


Xp(z) 

(z)  +  T 

z~'  Xs(2) 

(2-114) 

xp(z) 

■  z’’  xs 

(z)  +  T 

z'1  i\U) 

(2-115) 

xs(z) 

=  (l-o) 

xp(z)  + 

“  xm(z) 

(2-116) 

\(z) 

=  X  (z) 
p 

+  (S/T) 

xm(z)  -  (B/T)  Xp(z) 

(2-117) 

x-s(z) 

-1  •• 
=  Z  X 

S(z)  + 

(y/T2)  xm(z)  -  (y/T2)  xp(z)  . 

(2-118) 

In  the  alpha-beta-gamma  tracker,  there  is  one  input  to  the  system, 

x  (k),  and  five  output  signals  x  (k),  x  (k),  x  (k),  x  (k),  and  x’  (k). 
m  p  p  s  ^  s 

Therefore,  there  must  exist  five  separate  transfer  functions  relating 

x  (z)  to  x  (z),  x  (z),  x  (z),  x  (z),  and  x’  .  These  transfer  functions 
ni  p  p  s  s  s 

rnay  be  obtained  by  dividing  each  side  of  equations  (2-114)  through  (2-118) 
by  xm(z)  to  yield 


X_(z)  ,  X  (z)  ,  X  (z) 


(2-119) 


H2(z)  = 
H3(z)  = 


i(z) 

05T 


X$(z) 

y77 


^  *  T  z-'  UU 


x  [zj 
m  ' 


V*7 


('-»>  xp(z)  .  xm(z) 


(2-120) 


(2-121 ) 
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X  (z)  X  (z) 
H4(z>  =  yiT  ■ 


(6/T)  xm(z) 

xm(z) 

m 


(0/T)x  (z) 

X  (z) 

m  '  • 


(2-122) 


H5(z) 


(z) 


X  (z) 
m 


!  x‘s(z)  (y/T2)  xm(z)  (y/T2)  xp(z) 

X  (z)  +  X  ( z )  x_(z) 

m  m  m 

(2-123) 


After  incorporating  the  transfer  function  relationships 
x  (z)  x  (z)  x  (z) 

Vz>  =  ■  Vz>  =  ifazj  •  k3(2>  =  yil  • 


x  (z)  X*  (z) 

h4<2>  ■  ysr  •  h5(z)  ■  ^ 


and  rearranging  terms,  the  set  of  five  equations  becomes 


H-|  (z )  -  z'1  H3(z)  -  T  z’1  H4(z)  =  0 


H2(z)  -  z'1  H4(z)  -  T  z"1  H5(z)  =  0 


(1-a)  H1 (z)  -  H3(z)  =  -a 

(8/T)  H^Z)  -  H2(z)  +  H4(z)  =  &/T 
(y/T2)  H1  (z)  +  o-z"1)  H5(z)  =  y/T2 


(2-124) 

(2-125) 

(2-126) 

(2-127) 

(2-128) 


The  transfer  functions,  H..(z),  now  have  the  rational  form  N^(z)/D(z). 
The  general  determinant  D(z)  is  given  by 
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1, 

0, 

-z-\ 

-ir\ 

0 

0, 

1, 

0, 

-z-\ 

-TZ'1 

(1-a) 

0, 

-1, 

0, 

0 

(e/T), 

-1, 

0, 

1, 

0 

(y/T2), 

0, 

0, 

0, 

n-z_1) 

(2-129) 


which  upon  evaluation  yields  [12] 

D(z)  =  (-z)-3  [z3+(«+B-3)z2  +  (3+y-2a-g)z  +  (a-1)]  .  (2-130) 

By  using  the  theory  of  determinants  [12]  and  performing  considerable 
algebra,  the  transfer  functions  for  the  smoothed  output  variables  are 
[12], 


x  (z) 

H,(z)  =  — 


x  m 


o  o 

aZ°  +  (g-2a)z  +  (a+B-y)z _ 

Z3  +  (a+g-3  )Z2  +  (3+y-2a-S  )z  +  (a-1) 


(2-131  ) 


x  (z)  1 


8Z3  +  (y-2g)z2  4-  (b-y) 


'  z3  +  (a+B-3)z2  +  (3+Y~2a-B)z  +  (a-1  ) 


(2-132) 


x’(z)  v 

H  (z)  =  — =  X 


z3  -  2z2  ±  z _ 

xm^  T2  z3  +  (a+£5-3)z2  +  (3+y-2a-g)z  +  (a-1) 


(2-133) 


These  transfer  functions  have  the  same  denominator  polynomial  and 
therefore  have  the  same  set  of  poles.  Observing  the  transfer  function 
solutions  for  the  alpha-beta  filter  outputs  x$(z)  (2-56)  and  x$(z)  (2-57), 
it  is  interesting  co  note  that  the  solutions  (2-131)  through  (2-132)  do 
not  reduce  tc  those  of  the  simpler  alpha-beta  filter  solutions  for  y  =  0. 
This  can  also  be  seen  from  the  original  difference  equation  statements  of 
both  systems.  The  alpha -beta -gamma  filter  can  be  reduced  to  behave  as 
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the  alpha-beta  filter  in  practice  by  letting  y  =  0  and  requiring  the 
boundary  condition  x*s(0)  =  0.  The  point  to  be  made  here  is  that  the 
alpha-beta-gamma  filter  is  not  merely  a  simple  extension  of  the  alpha- 
beta  filter. 


Stability  of  the  Alpha-Beta- 


Gamma  Filter 


A  stable  tracking  system  requires  that  a  bounded  filter  input  pro¬ 
duces  a  bounded  filter  output.  In  the  Z  plane,  this  strictly  means  that 
the  poles  of  the  characteristic  equation  lie  within  the  unit  circle. 
There  are  many  tests  for  stability  but  the  Jury  stability  criterion 
[10]  will  be  applied  here  as  it  was  for  the  alpha-beta  filter.  First, 


the  denominator  of  the  transfer  function  is  written  as 


F(z )  =  (1 -a )  +  (3  -2a-g+y)z  +  (a+B-3)z2  +  {1 )z3  (2-134) 

where 

F(z)  »  a„  <-  a,z’  *  a/  +  a3z3 
From  (2-134)  the  Jury  array  obtained  is 


The  requirement  that  F(z)  has  all  of  its  roots  within  the  unit 
circle  is  satisfied  by  the  conditions: 

Condi tion  1 :  F(1 )  >  0 

Condition  2:  (-1)3  F(-l)  >  0 

Condition  3:  |a  1  <  a^ 

Condi tion  4 :  ' b^ j  >  j  b^ ! 


u*  j.a" 
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where 

b 

Condition  1: 

Condition  2: 


Condition  3: 
where , 


3  a  a  , 
0  m-k 


am  3. 

m  k 


m=3 


Z  3  +  (a+B-3)z3  +  (3+y-2a-8)z  -l  (a-1) 


>  0 


z=l 


(l)3  +  (a+8-3)(l)2  +  (3+y-2j-b) 0 )  +  (a-1  )  >  0 


l+a+8-3+3+a-2a-8+a-2>0 


a  >  0 


(2-136) 


Z3  +  (a  +  8-3)z3  +  (3+y-2a-8)z  +  (a-1)  I  >  0 

|z=-l 

M)3  +  (a+6-3)(-l)2  +  (3+y-2a-p)  (-1  )  +  (a-1)  >  0 

-1  +a+6-3-3-y+2a+8+a-l  >0 
4a  +  2f3-y-8>0 
4a  +  2f3  -y  >  8 

I  +  f  '  i  >  1  •  (2-137) 


<  a 


2' 


| a-1  |  <  1 

o  •  ,  •  : 


(2-13 


o ! 


AD-A110  863 
UNCLASSIFIED 


8E0R8IA  INST  OF  TECH  ATLANTA  CNflNEEftlN#  EXPERIMENT  ~ rrr  e/n  • 

"Bsawssar”1" 
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Condition  4:  j  bg|  >  | | , 
where, 

a0  a3 

bp  =  =  (a-1 )  -1  =  a^  -2a 

a3  a0 

ao  ai 

b2  =  =  (a-l)(a+B-3)  -  (3-2a-B+y) (1 ) 

a3  a2 

2 

b^  =  a  +  aB  -  2a  +  y  , 

therefore, 

2  .  2 
|a  -2a |  >  ja  +aB-2a+y| 

aB  >  y  ,  (2-139) 

since  0  <  a  <  2  , 

B  >  0  .  (2-140) 

The  relations  just  presented  are  necessary  and  sufficient  for  the  poles 
of  the  characteristic  equation  to  lie  within  the  unit  circle.  These 
relations  are  very  useful  when  manipulating  a,  B,  and  y  in  designing  a 
filter  to  satisfy  a  particular  set  of  specifications. 


Pole  Analysis  for  the  Alpha- 
Beta-fiamma  tracking  Filter 
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When  one  of  the  equations  (2-116)  through  (2-118)  are  solved  to 
obtain  one  of  the  outputs  xout(z)»  the  output  will  be  of  the  rational 
form  [9] 


xout 


(z) 


n  (z-w  ) 
i=l  1 


2  m. 

n  (z-p.) 

i=l  1 


(2-141) 


where  rtK  is  the  multiplicity  of  the  ith  pole,  wi  is  the  ith  zero  of  the 
system  and  pi  is  the  ith  pole  of  the  system.  The  partial  fraction  ex¬ 
pansion  and  inversion  will  yield  the  general  form  [9] 

ml 

*out<nT>  =  A0^pi,wi^  +  f.  ^(k,  Pi»  wi )  (P] ) 

IX™  l 


^  I 

+  s  ^  g2k(k,  pr  w.)(p2) 


+ 


mN 

2  k 

+  ^  9N2(k,  p.,  w.Mp^T  .  (2-142) 

All  the  terms,  except  Aq  decay  as  the  power  of  a  given  pole  or  root  of 

k 

the  denominator  (p.)  where  |p.|  <  1.  The  g  coefficients  may  contain 
powers  of  k  which  offset  this  decay,  but  the  structure  of  (2-142)  pro¬ 
vides  an  approach  to  alter  the  response  of  the  filter  by  choosing  the 
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poles  and  then  obtaining  the  required  a,  8,  and  y.  Thus  given  the  poles 
of  the  characteristic  equation  and  defining 


Dq(z)  =  Z3  +  (a+8~3)z2  +  (3+y-2a-8)z  +  (a-1 ) 


=  (z-p1){z-p2)(z-p3)  , 

(2-143) 

!t  can  be  shown  that  [12], 

P  =  a  +  8  -  3  =  -(p]  +  p2  +  p3) 

(2-144) 

Q  i  3  +  y  -  Za  ■  6  ■  p,p2  +  p2p3  *  P]p3 

(2-145) 

R  =  a  -  1  =  “P1P2P3  ’ 

(2-146) 

Conversely  the  poles  may  be  derived  by  knowing  a. 

8,  y  and  solving 

the  cubic  equation  DQ(z)  =  0.  If 

Dq(z)  =  z3  +  Pz2  +  Qz  +  R  , 

(2-147) 

then  from  [13], 

a  e  \  (3Q  -  P2) 

(2-148) 

be—  (2P3  -  9PQ  +  27R) 

(2-149) 

r  2  3"i  1/2 

G  5  [r 4  It] 

(2-150) 

(2-151 

"[*"]  '/3 

(2-152) 

=  A  +  B 

(2-153) 
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*2.3  ■-  [¥]±t*  [¥] 
and  the  poles  are, 

P]  =  E1  -  P/3 

p2,3  =  E2,3  '  P/3  • 

Alp ha -Beta -Gamma  Filter  Noise  Response 

Singer  and  Behnke  [14]  have  shown  that  an  alpha- beta  filter  has 
a  noise  response  which  is  less  than  20%  worse  than  that  of  a  Kalman  fil¬ 
ter  for  smoothed  estimates  of  position  and  velocity.  The  second  order 
comparison  suggests  that  the  alpha -beta-gamma  filter  will  perform  in  a 
similar  manner  when  compared  with  a  third  order  Kalman  filter. 

Once  again  the  variance  reduction  ratios  K  (0),  K*(0),  and  K - - (0) 

XX  X 

are  introduced  to  evaluate  the  noise  throughput  of  the  tracking  system. 
Simpson  [4]  defines  the  noise  performance  measures  as 

K  (0)  =  z  g  2(k)  (2-157) 

x  k=0  x 

K* (0)  *  K  z  g.2(k)  (2-158) 

x  T  k-0 

Kj.(0)-1t  z  g..2(k)  .  (2-159) 

x  r  k=0  x 

Wilcox  [12]  provides  an  algorithmic  approach  to  the  calculation  of  the 

noise  performance  of  the  alpha-beta-gamma  filter,  assuming  a  zero  mean 

2 

Gaussian  measurement  noise  with  variance  =  . 


(2-154) 

(2-155) 

(2-156) 
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For  the  calculation  of  Kx(0),  let 

bp  =  a,  b.j  =  0-2a,  bg  =  a+Y“0»  ^3  =  0 
30=1,  a  1  =  a+8- 3 ,  &2  =  3+y-2a-6*  3 3  =  a-1 

For  the  calculation  of  10(0),  let  the  a^s  remain  the  same  and 

bQ  =  (g/T),  b]  =  b2  =  b3  =  0  . 

For  the  calculation  of  K*«(0),  let  the  a^ 's  remain  the  same  and 
bQ  =  b2  =  (y/T2),  b1  =  (-2y/T2),  b3  =  0  . 

Now  defining 

Kx(0)  =  K-|  (0) ,  K*  (0)  =  K2(0),  K..(0)  =  K3(0), 

then  the  noise  response  for  each  estimate  is  calculated  by  substituting 
the  appropriate  values  of  a^  and  bi  into 


^(0)  = 


a0B0^O  ”  a0Bl^l  +  a0B2^2  "  B3^3 
a0^aQ2'a32^0  '  ^a0al"a2a3^1  +  ^a0a2'ala3^2  , 


(2-160) 


where 

B0  •  b02  ♦  bl2  +  b22  *  b32 
B]  *  2^b0bl  +  blb2  +  b2b3^ 
B2  *  2<b0b2  *  blb3> 
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E2  =  al  +  a3 

Q0  =  a0El  "  a3E2 
Q]  =  aQa1  -  a2a3 

Q2  =  a-|E2  -  a2E1 

^3  =  (aTa3)(E2  ~E1  ^  +  a0  ^a0E2‘a3El^ 


It  can  be  shown  that  as  a,  8.  and  y  approach  zero  the  poles  of  the 
characteristic  equation  approach  one,  which  results  in  higher  noise  sup¬ 
pression  since  the  bandwidth  decreases  and  the  effective  data  smoothing 
capability  of  the  filter  increases. 


Performance  Measures  for  the 
Alpha-Beta-Gamma  Filter 

Simpson  [4]  and  Neal  [15]  define  two  transient  performance  mea¬ 
sures  for  the  transient  analysis  of  the  alpha-beta-gamma  filter.  The 
first  transient  performance  measure  is  the  unit-increment  ramp  response 
used  earlier  in  the  analysis  of  the  alpha-beta  tracker, 


2 

z  [(unit-increment  ramp)  -  (position  ramp  response)] 
n=0 

(2-161) 


2 

x 


00 

Z  [(velocity  of  unit-increment  ramp) 
n=0 

2 

-  (velocity  ramp  response)] 


(2-162) 


1 

I 

I 

1 

;1 

\ 

j 
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OO 

D**  *  i  [(acceleration  of  unit  ramp  response  (=0)) 
x  n=0 

2 

-  (acceleration  ramp  response)]  ,  (2-463) 


thus, 

D2  =  E  [n  -  z  g  (j)(n-j)]2  (2-164) 

x  n=0  j=0  x 


D;2  =  E  [1  -  I  g-(j)(n-j)]2  (2-165) 

x  n=0  j=0  x 


D.-2=  L  [-  E  9**(j)(n-j)]2  .  (2-166) 

x  n=0  j=0  x 

2 

Similarly  the  transient  performance  measure  (A  )  due  to  a  unit  step  in¬ 
put  of  acceleration  is  defined,  in  the  same  manner,  as  the  sum  of  the 
squares  of  the  errors  arising  from  this  input.  Therefore, 


A  2  =  E  d--\  l  gx(j)(n-j)2]2  (2-167) 

x  n=0  c  c  j=0  x 


A-2  =4  £  [n  -  ]■  r  g-(j)(n-j)232  (2-168) 

x  T  n=0  c  ns0  x 

A-*2  =  E  [1-1  £  g *  *  ( j ) (n- j  )2]2  .  (2-169) 

x  T4  n=0  c  j*0  x 

By  numerical  evaluation  of  these  performance  measures,  Simpson 
[4]  derived  an  optimal  relation  between  a,  8,  and  y  such  that  for  a 
given  noise  smoothing,  the  optimal  transient  response  is  achieved,  and 
vice  versa. 

28  -  a(a  +  8  +  j  y)  =  0 


(2-170) 
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In  a  later  paper,  Neal  [15]  verified  Simpson's  relation  analyt¬ 
ically  using  a  new  approach  [1]  to  linear  estimation  theory  and  also 
derived  an  additional  relation  between  a,  6,  y. 

B2  =  2ay  (2-171) 

The  preceding  relations  (2-170)  and  (2-171)  provide  a  method  of  opti¬ 
mally  choosing  two  of  the  parameters  of  the  alpha-beta-gamma  tracking 
filter  thus  simplifying  the  task  of  optimizing  the  tracker  in  a  given 
tracking  environment. 

The  recursive  nature  of  the  alpha-beta-gamma  filter  makes  it  a 
prime  candidate  for  implementation  on  a  digital  computer.  A  flowchart 
of  the  alpha-beta-gamma  filter  is  given  in  Figure  2.5. 


III.  THE  TRI-STATE  ADAPTIVE  FILTER 
Introduction 

Three  different  tracking  filters  have  been  discussed  in  the  pre¬ 
vious  chapter  for  implementation  in  a  track-while-scan  radar  system. 
These  trackers  are  capable  of  providing  optimum  estimates  of  an  air¬ 
craft's  position  provided  the  dynamical  models  on  which  the  tracking  al¬ 
gorithm  are  based,  are  accurate  assumptions  of  the  aircraft's  true  tra¬ 
jectory.  Unfortunately,  due  to  the  presence  of  radar  measurement  noise 
and  wind  turbulence,  no  single  dynamical  model  for  the  aircraft’s  tra¬ 
jectory  will  be  an  accurate  model  throughout  the  period  that  the  air¬ 
craft  is  under  track. 


Tracking  Filter  Limitations 

Although  the  alpha  filter  is  the  simplest  tracker  that  has  been 
discussed,  it  provides  the  greatest  noise  smoothing  performance.  It  is 
based  upon  the  assumption  that  the  aircraft's  position  is  not  changing 
with  respect  to  time.  The  constant  position  model  is  applicable  to  the 
azimuth  or  elevation  tracking  of  an  aircraft  whose  trajectory  is  radial 
to  the  radar.  The  alpha  filter  is  capable  of  tracking  a  constant  posi¬ 
tion  target  with  less  error  than  either  the  alpha-beta  filter  or  the 
alpha- beta -gamma  filter,  but  because  the  alpha  filter  assumes  a  zero 
target  velocity,  its  ability  to  track  a  target  undergoing  a  change  in 
position  is  impaired. 
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The  alpha-beta  tracker  is  a  second  order  tracking  filter  which  is 
designed  to  minimize  the  mean-square  error  in  filtered  position  and  ve¬ 
locity,  assuming  no  velocity  change  between  data  samples.  Thus  the  alpha 
beta  filter's  capability  to  track  severely  maneuvering  (accelerating) 
targets  is  severely  limited. 

If  the  target  is  assumed  to  be  maneuvering,  its  trajectory  may  be 
modeled  by  a  constant  acceleration  model  with  time  varying  position  and 
and  velocity.  This  more  generalized  model  of  the  target's  trajectory 
leads  to  the  third  order  alpha-beta-gamma  tracker  which  is  capable  of 
providing  good  estimates  of  the  target's  position  and  velocity  through¬ 
out  a  maneuver.  If,  however,  the  target  is  not  maneuvering,  the  alpha- 
beta-gamma  filter  suffers  a  significant  degradation  in  noise  smoothing 
performance  compared  to  the  simpler  alpha-beta  tracker  (for  constant 
velocity  traj ■ /'tory)  or  alpha  tracker  (for  constant  position  trajectory). 

Selection  of  the  Appropriate  Filter 

Since  the  state  of  the  target's  trajectory  may  change  as  a  func¬ 
tion  of  time  and  range,  no  single  filter  will  provide  accurate  estimates 
of  the  target's  position  and  velocity  throughout  the  tracking  interval. 
However,  if  a  bank  of  three  filters  consisting  of  the  alpha  filter,  the 
alpha-beta  filter,  and  the  alpha -beta -gamma  filter  were  operated  in  par¬ 
allel  with  one  another,  accurate  estimates  of  the  target's  position  and 
velocity  would  be  available  regardless  of  the  changes  in  the  target's 
trajectory.  Some  type  of  adaptivity  should  be  built  into  the  filter  bank 
so  that  if  no  position  change  is  detected  then  the  estimate  from  the 
alpha  tracker,  should  be  selected.  If  the  target  is  ascertained  to  have 
a  constant  velocity,  then  the  estimate  from  the  alpha-beta  tracking 
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filter  should  be  selected.  And  if  the  target  is  undergoing  a  change  in 
velocity,  then  the  estimate  from  the  alpha -beta -gamma  tracker,  should 
be  selected. 

Detection  of  a  change  in  the  target's  trajectory  and  the  selection 
of  the  appropriate  filter  output  may  be  accomplished  by  monitoring  the 
variance  in  the  smoothed  position  estimates  of  both  the  alpha  filter  and 
the  al pha-beta  f il ter. 

Alpha  Filter  Output  Variance  as  a 
Basis  for  Filter  Selection 

If  a  constant  position  measurement  of  the  form 

xm(k)  «  X  +  n(k)  (3-1) 

where  , 

x  (k)  -  radar  position  measurement  at  time  kT 
m 

X  =  actual  position  at  time  kT  (constant) 

n(k)  =  zero  mean  Gaussian  noise  with  variance  =  o  > 

n 

is  the  input  to  an  alpha  tracker,  then  from  equation  (2-34)  the  variance 
in  the  smoothed  position  estimate  is  given  by 

varlx  (k)}  =  o  2  .  (3-2) 

s  z-a  n 

2 

The  input  noise  variance  on  is  usually  available,  and  if  not,  it  can  be 
easily  calculated.  The  variance  in  the  smoothed  position  estimates  of 
the  alpha  filter  can  be  approximated  directly  from  the  filter's  output 
by 
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var{xs(k)}  -  |- 


L  2 

r  [x  (1)  -  E{x  (k)}] 

1  =  0  s  s 


(3-3) 


where, 

*s(i)  =  ith  estimate  of  smoothed  position 

L  =  the  integral  number  of  samples  over  which  the 
variance  in  smoothed  position  is  being  calculated 

!  L 

E{x  (k) }  =  r  2  x  (i )  . 

s  L  i=0  s 

As  long  as  constant  position  measurements  of  the  form  prescribed 
in  equation  (3-1),  are  input  to  the  alpha  tracker,  then  the  variance  in 
the  smoothed  position  estimates,  approximated  by  equation  (3-3),  will 
remain  less  than  or  equal  to  the  variance  calculated  via  equation  (3-2). 
If,  however,  radar  measurements  of  a  target  undergoing  a  change  in  posi¬ 
tion  are  input  to  the  alpha  tracker,  the  variance  in  the  smoothed  posi¬ 
tion  estimates,  approximated  by  equation  (3-3),  will  exceed  the  expected 
variance  calculated  via  equation  (3-2).  Thus  the  variance  in  smoothed 
position  obtained  from  equation  (3-2)  may  be  used  as  a  threshold  to 
determine  if  the  alpha  filter  is  capable  of  tracking  the  target.  If  the 
variance  in  the  smoothed  position  estimates,  approximated  by  equation 
(3-3),  is  less  than  the  variance  threshold  calculated  via  equation  (3-2), 
then  the  alpha  tracker's  smoothed  estimate  of  position  is  considered  the 
best  estimate  of  the  target's  position  and  the  velocity  of  the  target  is 
assumed  to  be  zero.  If  the  variance  of  the  alpha  tracker's  smoothed  po¬ 
sition  estimates,  approximated  by  equation  (3-3),  exceeds  the  variance 
threshold  obtained  from  equation  (3-2),  then  the  alpha  tracker's  smoothed 
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estimate  of  position  is  rejected  and  the  alpha-beta  filter's  smoothed  es¬ 
timates  of  position  and  velocity  are  examined  to  determine  if  they  yield 
a  more  accurate  representation  of  the  target's  trajectory.  When  the 
variance  of  the  alpha  tracker’s  smoothed  output  (3-3)  drops  below  the 
variance  threshold  (3-2),  then  the  alpha  filter's  smoothed  position  es¬ 
timate  is  once  again  considered  a  valid  estimate  of  the  target's  position. 


Alpha-Beta  Filter  Output  Variance 
as  a  Basis  for  Filter  Selection 

In  a  similar  manner,  if  the  radar  measurement  of  a  constant  velo¬ 
city  target  is  input  to  an  alpha-beta  filter  it  has  the  form  of 

xm(k)  =  x(k)  +  vkT  +  n(k)  ,  (3-4) 


where, 

xm(k)  *  radar  position  measurement  at  time  kT 
x(k)  *  actual  position  at  time  kT 
v  *  velocity  of  target  (constant) 

T  *  sampling  interval 

Assuming  an  input  as  prescribed  in  (3-4),  the  variance  expected  in  the 
smoothed  position  estimate  is 


var{x$(k)}  = 


2a2  +  s(2-3a)l  2 

a(4-S-2a)  an 


(3-5) 


Equation  (3-5)  is  of  course,  based  upon  the  previous  relation  given  by 
equation  (2-91). 

As  with  the  alpha  filter,  the  variance  in  the  smoothed  position 
estimates  of  the  alpha-beta  filter  can  be  approximated  from  the  filter's 
smoothed  output  by  equation  (3-3). 
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As  long  as  the  target  maintains  a  constant-velocity  straight-line 
trajectory,  the  variance  in  the  alpha-beta  filter's  smoothed  estimates 
of  position  (3-3)  will  remain  less  than  or  equal  to  the  expected  vari¬ 
ance  in  smoothed  position  calculated  via  equation  (3-5).  Since  the 
alpha-filter's  estimate  of  position  has  been  previously  examined  and 
found  to  be  inaccurate  for  a  constant-velocity  target  trajectory,  it  is 
assumed  that  the  alpha-beta  filter's  smoothed  estimates  of  position  and 
velocity  best  represent  the  true  trajectory  of  the  target. 

If  the  target  undergoes  an  acceleration,  the  variance  in  the  alpha- 
beta  filter's  smoothed  position  estimates,  approximated  by  equation  (3-3), 
will  exceed  the  variance  threshold  calculated  in  equation  (3-6).  Under 
this  condition,  the  alpha-beta  filter's  smoothed  estimates  of  position 
and  velocity  are  deemed  inaccurate  and  tracking  control  is  transferred 
to  the  higher  order  alpha-beta-gamma  filter.  Estimates  of  smoothed  posi¬ 
tion  and  smoothed  velocity  from  the  alpha-beta-gamma  filter  are  used  only 
as  long  as  the  variances  in  the  smoothed  position  estimates  of  the  alpha 
filter  or  alpha-beta  filter  exceed  their  variance  thresholds,  calcu¬ 
lated  via  equations  (3-2)  and  (3-5),  respectively.  If  the  variances  in 
the  smoothed  position  estimates  of  both  the  alpha  filter  and  the  alpha- 
beta  filter  drop  below  their  respective  variance  thresholds  during  the 
same  sample  interval,  then  track  control  is  transferred  to  the  alpha  fil¬ 
ter  since  it  provides  a  greater  degree  of  noise  smoothing  than  the  other 
two  filters.  A  block  diagram  depicting  the  tracking  system  discussed  in 
the  previous  paragraphs  is  shown  in  Figure  3.1. 


Figure  3.1.  Block  Diagram  of  the  Tri-State  Adaptive  Tracking  Filter 
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Implementation  of  the  Tri-State  Adaptive 
Tracking  Filter 

Examination  of  Figure  3.1  shows  that  the  alpha,  alpha-beta,  and 
alpha-beta-gamma  tracking  filters  operate  in  parallel.  Each  filter  gen¬ 
erates  its  own  estimates  of  smoothed  position  and  smoothed  velocity 
based  upon  its  implicit  model  of  the  target's  trajectory.  The  smoothed 
position  estimates  from  the  alpha  filter  and  the  alpha-beta  filter  are 
stored  in  separate  L-length  shift  registers.  From  the  stored  values  of 
the  alpha  filter's  smoothed  position  estimates,  the  variance  in  the  alpha 
filter's  smoothed  position  estimates  is  calculated  via  equation  (3-3). 

In  a  similar  manner  the  variance  in  the  alpha-beta  filter's  smoothed 
position  is  calculated  using  the  stored  values  of  the  alpha-beta  filter's 
smoothed  position  estimates.  The  calculated  variances  in  the  alpha  fil¬ 
ter  and  alpha-beta  filter  smoothed  position  estimates  are  then  compared 
to  their  respective  variance  thresholds  (calculated  via  (3-2)  and  (3-5)). 
The  appropriate  filter  output  is  then  selected  based  upon  the  criteria 
previously  discussed.  At  the  next  sample  interval,  the  new  smoothed 
position  estimates  are  shifted  into  the  shift  registers  and  the  variance 
in  smoothed  position  is  recalculated  for  each  of  the  two  filters. 

A  flowchart  depicting  the  operation  of  the  tri-state  adaptive 
tracking  algorithm  is  presented  in  Figure  3.2. 

The  implementation  of  the  tri-state  adaptive  tracking  filter,  as  a 
FORTRAN  IV  subroutine,  is  given  in  the  Appendix.  The  subroutine  per¬ 
forms  azimuth  tracking  on  an  incoming  aircraft.  The  inputs  to  the  sub¬ 
routine  are  the  sampling  interval  T,  the  azimuth  noise  variance,  the 
range  of  the  aircraft  at  each  sample  interval,  and  the  radar's  measure¬ 
ment  of  lateral  position  in  feet. 


IV.  SYSTEM  SIMULATION 


Introduction 

This  study  investigates  the  application  of  a  tri-state  adaptive 
filter  in  an  automatic  landing  system  for  aircraft.  The  automatic  land¬ 
ing  system  under  consideration  is  the  Marine  Air  Traffic  Control  and 
Landing  System  (MATCALS).  The  MATCALS  utilizes  an  AN/TPN-22  track-while- 
scan,  electronically  steerable  microwave  radar  to  produce  raw  digitized 
measurements  of  the  aircraft's  vertical  and  lateral  positions.  This  po¬ 
sitional  information  is  filtered  and  processed  by  a  groundbased  control¬ 
ler  which  calculates  the  appropriate  bank  and  pitch  commands  to  correct 
the  aircraft's  trajectory.  The  commands  are  then  transmitted  via  a  com¬ 
munication  link  to  the  aircraft's  autopilot  to  automatically  land  the 
aircraft.  This  configuration  is  presented  in  Figure  4.1.  A  more  de¬ 
tailed  discussion  of  the  MATCALS  is  available  in  [16]. 

FORTRAN  IV  programs  have  been  written  to  simulate  the  F4J  and  A7E 
aircraft  dynamics,  the  F40  and  A7E  autopilots,  the  landing  system  con¬ 
troller,  and  the  AN/TPN-22  radar.  The  FORTRAN  IV  program  listings  are 
given  in  [16].  Each  of  the  aircraft  simulations  consists  of  two  parts  - 
the  lateral  control  system  simulation  and  the  longitudinal  control  system 
simulation.  The  two  control  systems  are  uncoupled  .  are  structurally 
identical.  The  lateral  control  system  of  the  F4J  ,  . i  be  used  to  study 
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Figure  4.1.  Block  Diagram  of  MATCALS  Automatic  Landing  System. 
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the  performance  of  the  tri-state  adaptive  filter,  since  stability  is 
more  difficult  to  obtain  in  this  system  than  in  the  longitudinal  control 
system  [17]. 


The  F4J  Control  System 

The  lateral  control  system  of  the  F4J  automatic  landing  system  is 
modeled  as  a  sampled  data  control  system.  The  aircraft's  lateral  dynam¬ 
ics  and  autopilot  generate  the  analog  signal  y(t)  (lateral  position). 

The  AN/TPN-22  radar  is  modeled  as  an  A/D  converter,  which  converts  the 
continuous  process  y(t)  into  discrete  measurements  of  position  every  T 
seconds,  where  T,  the  sampling  period,  is  selected  to  be  0.1  seconds. 

The  groundbased  controller  is  modeled  as  a  digital  transfer  function 
which  generates  discrete  bank  commands,  <t>(k) ,  and  a  zero  order  hold 
which  restructures  the  discrete  bank  commands  into  the  continuous  sig¬ 
nals  <(>(t).  A  block  diagram  of  the  lateral  control  system  for  the  F4J 
aircraft  is  given  in  Figure  4.2. 

The  F4J  Lateral  Guidance  System 

The  F4J  aircraft  lateral  dynamics  are  described  by  a  set  of  sixth 
order  linear  differential  equations.  The  autopilot  dynamics  are  given  by 
a  set  of  third  order  differential  equations.  Thus,  the  lateral  control 
system  is  ninth  order.  In  addition  the  lateral  control  system  includes 
three  nonlinearities.  The  system  equations  are  of  the  form 


x(t)  =  Ax(t)  +  Bu(t)  +  Ef(t) 

(4-1) 

£(t)  =  Cx(t)  +  Du(t) 

(4-2) 
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Figure  4.2.  Block  Diagram  of  the  F4J  Lateral  Control  System. 
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where  the  nonlinearities  are  simulated  by  the  term  Ef .  The  simulation  of 
the  F4J  aircraft's  lateral  dynamics  and  the  dynamics  of  the  autopilot  is 
discussed  in  detail  in  Reference  [16]. 

The  SPN-42  Digital  Controller 

The  groundbased  control  unit  for  the  F4J  lateral  control  system  is 
the  SPN-42  digital  controller.  The  SPN-42  is  described  as  a  PID  control¬ 
ler.  In  this  notation,  P  indicates  proportional,  I  indicates  integral, 
and  D  indicates  derivative.  The  SPN-42  digital  controller  is  composed  of 
an  a-e  tracking  filter,  a  differentiator,  an  integrator,  four  a  filters, 
and  a  floating  limiter.  A  block  diagram  of  the  SPN-42  lateral  control 
system  digital  controller  is  presented  in  Figure  4.3,  with  all  nonlinear¬ 
ities  omitted.  A  detailed  description  of  the  SPN-42  digital  controller 
is  available  in  [17]  . 

AN/TPN-22  Radar 

The  AN/TPN-22  phased  array  radar  is  utilized  in  the  lateral  control 
system  of  the  F4J  aircraft  to  measure  the  lateral  position  of  the  air¬ 
craft.  In  the  measurement  process  a  significant  amount  of  noise  is  in¬ 
troduced  into  the  system  such  that  the  reported  lateral  position  of  the 
aircraft  is  degraded.  A  model  for  simulating  the  noise  associated  with 
the  AN/TPN-22  radar  is  the  ITT-Gilfillan  AN/TPN-22  radar  noise  model 
[18].  A  block  diagram  of  the  radar  noise  simulation  is  given  in  Fig¬ 
ure  4.4.  The  difference  equation  describing  the  AN/TPN-22  azimuth  radar 
noise  is  given  by  [17], 


A<t(n)  =  0.382  A*(n-1)  +  0.15  t$(n-2)  +  0.122  «>(n-3) 

+  0.0045  A  (n-4)  +  v(n) 


(4-3) 
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Figure  4.3.  SPN-42  Digital  Controller. 
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Figure  4.4.  AN/TPN-22  Radar  Error  Simulation. 
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and, 

measured  angle  (n)  =  A<p(n)  +  g(n)  +  target  angle  ,  (4-4) 

where, 

A$(n)  =  azimuth  noise  at  the  nth  sample 
g(n)  =  azimuth  granularity  error  sequence  with  standard 
deviation,  =  0.148  mr 

v(n)  =  white  noise  sequence  with  mean  =0.0  and  variance, 
c/  =  0.320  mr2 

The  aircraft's  normal  touchdown  point  on  the  runway  is  offset  from  the 
radar  center  by: 

Longitudinal  X  offset  =  762.9  feet 

Lateral  Y  offset  =  -178.1  feet 

The  target  angle  is  given  by: 

target  angle  =  tan"1  £  x'T'762’9  ]  '  (4-5) 

The  lateral  position  of  the  aircraft,  measured  by  the  radar,  is  given  by 
YR(n)  =  (X  +  762.9)  tan[measured  angle  (n)]  -  178.1  (4-6) 


V.  SIMULATION  RESULTS 


Introduction 

The  feasibility  of  using  the  tri-state  adaptive  tracking  filter  in 
the  MATCALS  may  be  studied  by  incorporating  the  tri-state  adaptive  filter 
and  the  alpha-beta  filter  into  the  F4J  lateral  control  system  simulation 
discussed  in  Chapter  IV.  The  response  of  the  F4J  lateral  control  system, 
using  the  alpha-beta  filter  to  estimate  the  aircraft's  position  and  ve¬ 
locity,  may  then  be  compared  to  the  control  system's  response  using  the 
tri-state  adaptive  filter  to  estimate  the  aircraft's  position  and  velo¬ 
city.  Block  diagrams  depicting  the  implementation  of  the  alpha-beta 
filter  and  the  tri-state  adaptive  filter  in  the  F4J  lateral  control  sys¬ 
tem  are  given  in  Figure  5-1  and  5-2,  respectively. 

For  ease  of  reference,  the  lateral  control  system  using  the  alpha- 
beta  filter  to  estimate  the  aircraft's  position  and  velocity  shall  be 
called  the  alpha-beta  control  system.  Likewise,  the  lateral  control 
system  using  the  tri-state  adaptive  filter  to  estimate  the  aircraft's 
position  and  velocity  shall  be  called  the  tri-state  adaptive  control 
system. 


Tracking  System  Comparison 

Filter  Parameters 

The  filter  coefficients  used  to  obtain  the  simulation  results  for 
the  alpha-beta  control  system  and  the  tri-state  adaptive  control  system 


are  given  in  Table  5-1.  The  filter  coefficients  for  the  alpha-beta 
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control  system  were  obtained  from  [17],  and  were  found  to  provide  the 
optimal  noise  smoothing  performance  for  a  given  maneuver  following  capa¬ 
bility.  For  this  reason,  the  alpha-beta  stage  of  the  tri-state  adaptive 
filter  is  identical  to  the  filter  used  in  the  alpha-beta  control  system. 
The  filter  coefficients  for  the  alpha  and  alpha-beta-gamma  stages  of  the 
tri-state  adaptive  filter  were  determined  experimentally  by  making  re¬ 
petitive  simulation  runs  and  selecting  the  coefficients  which  provided 
the  best  noise  smoothing  performance  (alpha-filter)  and  the  best  maneu¬ 
ver  following  capability  (alpha-beta-gartma  filter). 

Filter  Frequency  Response 

A  FORTRAN  program  from  [19]  was  used  to  calculate  the  amplitude 
and  phase  responses  of  the  alpha-beta  filter  and  the  three  filters  which 
make  up  the  tri-state  adaptive  filter.  These  frequency  responses  are 
given  in  Figure  5-3  through  Figure  5-6. 

Control  System  Time  Response  with 
a  Given  Initial  Condition 

The  initial  condition  time  responses  of  the  alpha-beta  control 
system  and  the  tri -state  adaptive  control  system  are  given  in  Figure 
5-7  and  Figure  5-8. 

The  time  responses  of  the  two  lateral  control  systems  were  ob¬ 
tained  from  simulation  runs  of  the  eighty  seconds  of  the  aircraft's 
flight  prior  to  touchdown.  This  corresponds  to  an  initial  range  of 
17,632  feet.  The  aircraft  is  initially  assumed  to  be  offset  laterally 
from  the  extended  centerline  of  the  runway  by  twenty  feet.  The  lateral 
velocity  of  the  aircraft  is  initially  assumed  to  be  zero  while  the 
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forward  velocity  of  the  aircraft  is  assumed  to  be  a  constant  220.39 
feet  per  second.  For  the  time  response  simulations,  the  wind  and  noise 
disturbances  were  eliminated. 

The  time  responses  of  the  alpha-beta  control  system  and  the  tri- 
state  adaptive  control  system  may  be  evaluated  by  comparing  the  tran¬ 
sient  response  of  each  system  and  the  ability  of  each  system  to  guide 
the  aircraft.  The  rise  time  and  the  time  to  peak  provide  a  measure  of 
the  transient  responses  of  the  lateral  control  systems.  The  rise  time, 

T  ,  is  defined  as  the  time  required  for  the  aircraft  to  reach  the  ex¬ 
tended  centerline  of  the  runway.  The  time  to  peak,  Tp,  is  defined  as 
the  time  required  for  the  aircraft  to  reach  its  first  peak.  An  idea  of 
how  well  the  lateral  control  systems  guide  the  aircraft  may  be  obtained 
by  examining  the  percent  overshoot  and  the  time  to  settle  for  each  sys¬ 
tem.  The  percent  overshoot  for  the  lateral  control  system  is  given  by 

X  . 

P.0.  =  -j—1  x  100%  (5-1) 

*IC 


where  Xpeaj<  is  the  magnitude  of  the  first  peak  position  of  the  aircraft, 
and  is  the  magnitude  of  the  initial  position  of  the  aircraft.  The 

time  to  settle  is  defined  as  the  amount  of  time  necessary  for  the  con¬ 
trol  systems  to  settle  the  aircraft  within  2%  of  its  initial  position  of 
twenty  feet.  These  time  response  performance  measures  for  the  alpha- 
beta  control  system  and  the  tri -state  adaptive  control  system  are  given 
in  Table  5-2. 
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TABLE  5-2 

TIME  RESPONSE  CHARACTERISTICS 


ALPHA-BETA 
FILTER  CONTROL 
SYSTEM 

TRI-STATE 
ADAPTIVE 
CONTROL  SYSTEM 

TIME  TO  RISE,  T 

r 

( seconds ) 

9.8 

10.8 

TIME  TO  PEAK,  T 

P 

(seconds) 

18.3 

19.6 

TIME  TO  SETTLE,  T 

s 

(seconds ) 

56 

56 

PERCENT  OF  OVERSHOOT,  P.O. 

(*) 

36.75 

32.3 
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The  results  given  in  Table  5-2  reveal  that  the  alpha-beta  control 
system  has  a  faster  response  than  the  tri-state  adaptive  control  system, 
but  the  tri-state  adaptive  control  system  does  not  overshoot  as  severely. 

Tri-State  Adaptive  Control 
System  Performance 

The  performance  of  the  tri -state  adaptive  tracking  filter  may  best 
be  evaluated  by  comparing  the  results  of  a  number  of  simulation  runs  of 
the  tri-state  adaptive  control  system  with  the  results  of  similar  runs 
of  the  alpha-beta  control  system.  The  characteristics  of  the  system 
responses  cannot  be  obtained  from  one  simulation.  To  obtain  an  accurate 
statistical  description,  it  is  necessary  to  make  many  simulation  runs 
and  to  statistically  average  the  results. 

Monte  Carlo  Simulation  Runs 

One  method  of  obtaining  accurate  statistical  descriptions  of  the 
two  lateral  control  system  responses  is  by  making  Monte  Carlo  runs.  Each 
of  the  Monte  Carlo  runs  that  are  to  be  presented  is  determined  from 
twenty  simulation  runs.  For  each  simulation  run,  the  aircraft  is  ini¬ 
tially  assumed  to  be  on  the  extended  centerline  of  the  runway  at  a 
starting  range  of  17,632  feet.  Prior  to  the  start  of  each  simulation 
run,  the  aircraft  is  assumed  to  have  zero  lateral  velocity,  and  a  con¬ 
stant  forward  velocity  of  220.39  feet  per  second.  The  forward  velocity 
of  220.39  feet  per  second  results  in  a  flight  time  of  eighty  seconds  to 
touchdown  point.  All  lateral  movements  of  the  aircraft  during  the  sim¬ 
ulation  runs  are  caused  by  wind  turbulence  or  radar  noise.  The  wind 
turbulence  and  radar  noise  for  each  of  the  twenty  simulation  runs  have 
the  same  statistical  parameters;  however ,  the  random  number  generators 
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in  the  wind  turbulence  and  radar  noise  sources  are  started  at  different 
values  thus  generating  different  number  sequences  for  each  simulation 
run.  It  should  be  noted  that  each  Monte  Carlo  run  uses  identical  sets 
of  random  number  sequences. 

For  the  purpose  of  analyzing  the  two  control  system  responses  as 
a  function  of  range,  each  simulation  run  is  divided  into  four  equal 
time  bins,  each  having  a  duration  of  twenty  seconds.  This  corresponds 
to  four  range  bins,  each  having  a  length  of  4408  feet. 

Statistics  are  calculated  to  describe  the  system  responses  over 
the  appropriate  range  bin  for  each  simulation  run.  The  statistics  for 
each  range  bin  are  then  averaged  over  the  twenty  simulation  runs  which 
make  up  the  Monte  Carlo  run.  The  results,  as  a  function  of  range,  for 
the  Monte  Carlo  runs  with  wind  turbulence  and  radar  noise  included  are 
given  in  Table  5-3.  The  statistics  used  to  evaluate  the  alpha-beta  con¬ 
trol  system  and  the  tri-state  adaptive  control  system  are  the  mean 
square  lateral  position  error  off  the  extended  centerline  of  the  runway 
and  the  variance  in  position  error  off  the  extended  centerline  of  the 
runway.  A  listing  of  percent  improvement  of  the  tri -state  adaptive 
control  system  over  the  alpha-beta  control  system  is  included  in  the 
Tables.  Percent  improvement  is  defined  as 


%  Improvement  = 


(W 

TTma+mt)/2]  x  1 00% 


(5-2) 


where  M^  and  M^.  are  the  mean  square  error  of  position  off  the  extended 
centerline  of  the  runway  for  the  alpha-beta  control  system  and  the  tri- 
state  adaptive  control  system,  respectively. 


RESULTS  OF  MONTE  CARLO  RUNS 
WITH  WIND  AND  RADAR  NOISE 


81 


Table  5-3  indicates  that  the  tri-state  adaptive  control  system 
does  not  perform  as  well  as  the  alpha-beta  control  system  when  radar 
noise  and  wind  turbulence  is  included  in  the  simulation  runs.  Closer 
scrutiny  of  Table  5-3  reveals  that  the  tri -state  adaptive  control  system 
performs  best  only  for  bin  1  which  corresponds  to  a  range  greater  than 
13,224  feet. 

A  possible  reason  for  the  inferior  performance  of  the  tri-state 
adaptive  control  system  may  lie  in  the  switching  process  between  the 
alpha  and  alpha-beta  filter  states  of  the  tri-state  adaptive  filter. 

As  pointed  out  in  Chapter  III,  the  alpha  filter  will  provide  excellent 
smoothing  of  position  measurements  so  long  as  the  lateral  velocity  of 
the  aircraft  is  equal  to  zero.  If  the  aircraft  acquires  a  lateral  velo¬ 
city,  then  the  alpha  filter  estimate  will  be  in  error.  In  order  to 
determine  if  the  aircraft  is  acquiring  a  lateral  velocity,  the  variance 
in  the  smoothed  position  estimates  of  the  alpha  filter  state  must  be 
calculated  as  given  in  equation  (3-3).  This  variance  in  smoothed  posi¬ 
tion  is  then  compared  to  the  variance  threshold  for  the  alpha  filter 
given  by  equation  (3-2).  If  the  variance  in  the  alpha  filter's  smoothed 
position  estimates  exceeds  the  threshold  then  the  aircraft  is  assumed 
to  have  acquired  a  lateral  velocity  and  the  output  of  the  alpha-beta 
filter  should  be  selected.  If,  however,  the  variance  in  the  alpha 
filter's  smoothed  position  estimates  does  not  increase  rapidly  enough, 
the  tri-state  adaptive  filter  will  output  the  smoothed  position  estimate 
of  the  alpha  filter,  even  though  the  aircraft  has  acquired  a  lateral 
velocity,  and  thus  be  in  error. 


TABLE  5-4 


TABLE  5-5 
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simulation  runs  with  identical  wind  and  radar  noise.  One  run  uses  the 
alpha-beta  control  system  to  guide  the  aircraft  wb'le  the  second  run 
uses  the  tri-state  adaptive  control  system.  Figure  5-9  provides  a  com¬ 
parison  of  the  performance  of  the  two  control  systems  under  the  identi¬ 
cal  effects  of  wind  and  radar  noise.  The  tri -state  control  system  is 
seen  to  perform  better  than  the  alpha-beta  control  system,  especially 
during  the  first  fifty  seconds  of  the  simulation  run. 

The  improvement  in  performance  that  the  tri -state  adaptive  control 
system  has  over  the  alpha-beta  control  system  is  due  to  the  tri-state 
adaptive  filter's  more  accurate  estimates  of  velocity.  Shown  in  Figure 
5-10  is  a  comparison  of  the  actual  aircraft  lateral  velocity  to  the 
alpha-beta  filter  estimates  of  velocity.  A  comparison  of  the  actual 
aircraft  lateral  velocity  to  the  tri-state  adaptive  filter  estimates  is 
shown  in  Figure  5-11.  A  tabular  comparison  of  the  mean  square  error 
in  the  two  filter's  velocity  estimates  is  given  in  Table  5-7.  The  tri- 
state  adaptive  filter's  ability  to  more  accurately  estimate  the  air¬ 
craft's  lateral  velocity  is  the  key  to  the  improvement  in  the  performance 
of  the  tri-state  adaptive  control  system  over  the  alpha-beta  control  sys¬ 
tem  presently  used  in  the  system. 
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Comparison  of  the  Tri -State  Adaptive  Control  System  Response  to  the  Alpha-Beta  Control 
System  Response  with  Identical  Wind  and  Radar  Noise. 
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VI.  CONCLUSIONS 


Three  digital  tracking  filters,  each  based  upon  a  different  air¬ 
craft  dynamical  model,  were  combined  to  form  the  tri-state  adaptive 
tracking  filter.  The  selection  of  the  appropriate  filter  output  was 
determined  by  the  variance  of  the  filters'  smoothed  position  estimates. 
The  tri-state  adaptive  filter  was  implemented  in  the  simulation  of  the 
F4J  lateral  control  system.  The  results  given  in  Chapter  V  suggest  that 
the  performance  of  the  F4J  lateral  control  system  may  be  improved 
through  the  use  of  a  tri-state  adaptive  tracking  filter.  Since  the  F4J 
longitudinal  control  system  is  structurally  identical  to  the  lateral 
control  system,  the  tri-state  adaptive  tracking  filter  may,  in  a  similar 
manner,  provide  an  improvement  in  the  performance  of  the  longitudinal 
control  system. 

The  overall  performance  of  the  tri-state  adaptive  tracking  filter 
may  be  enhanced  by  selecting  the  parameters  of  each  of  the  three  compo¬ 
nent  filters  in  such  a  manner  as  to  achieve  a  more  complementary  filter 
response.  Another  modification  which  might  improve  the  performance  of 
the  tri-state  adaptive  filter  is  the  adjustment  of  the  variance  thres¬ 
holds  of  the  alpha  and  alpha-beta  filters.  As  was  shown  by  the  results 
of  the  F4J  lateral  control  system  simulation,  the  frequency  response  of 
the  tri -state  adaptive  filter  may  be  altered  by  the  selection  of  the  ap¬ 
propriate  variance  thresholds. 
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SUBROUTINE  FILBNK  (RADR) 

COMMON/ RANG/ RANGE 
COMMON/ FI LCON/SMP.SMV 

DIMENSION  SP1 (8) ,  SP2{8),  D1F1(8),  D1F2(8) 

DATA  SPl/8*0. 0/ , SP2/8*0. 0/.D1F1 /8*0. 0/ , D1 F2/8*0. 0/ 

CONVERT  INPUT  NOISE  VARIANCE,  VRP,  FROM  RADIANS  TO  FEET. 

VRP=4. 225E-07*(RANGE**2) 

SET  L,  THE  LENGTH  OF  THE  SHIFT  REGISTERS. 


L=5 

L2=L-1 

C 

C  TT= SAMPLE  INTERVALS  SEC. 

C  AT  IS  THE  POSITION  SMOOTHING  CONSTANT  OF  THE  ALPHA  FILTER. 

C  A2  IS  THE  POSITION  SMOOTHING  CONSTANT  OF  THE  ALPHA-BETA  FILTER. 

C  B2  IS  THE  VELOCITY  SMOOTHING  CONSTANT  OF  THE  ALPHA-BETA  FILTER. 

C  A3  IS  THE  POSITION  SMOOTHING  CONSTANT  OF  THE  ALPHA-BETA-GAMMA  FILTER. 

C  B3  IS  THE  VELOCITY  SMOOTHING  CONSTANT  OF  THE  ALPHA-BETA-GAMMA  FILTER. 

C  G  IS  THE  ACCELERATION  SMOOTHING  CONSTANT  OF  THE  ALPHA-BETA-GAMMA  FILTER. 
C 

C  RADR=RADAR  MEASUREMENT  OF  AZIMUTH  INPUT  TO  SUBROUTINE. 

C  SPI (L)=L-LENGTH  ARRAY  CONTAINING  THE  LAST  L-VALUES  OF  THE  ALPHA 
FILTER'S  SMOOTH  POSITION  ESTIMATES. 

C  SP2(L)-L  LENGTH  ARRAY  CONTAINING  THE  LAST  L-VALUES  OF  THE  ALPHA-BETA 
FILTER'S  SMOOTH  POSITION  ESTIMATES. 

C  SP3=ALPHA-BETA-GAMMA  FILTER'S  SMOOTH  POSITION  ESTIMATE. 

C  SV1 =ALPHA  FILTER'S  SMOOTH  VELOCITY  ESTIMATE 
C  SV2=ALPHA-BETA  FILTER'S  SMOOTH  VELOCITY  ESTIMATE. 

C  SV3=ALPHA-BETA-GAMMA  FILTERS  SMOOTH  VELOCITY  ESTIMATE. 

C  SA3=ALPHA- BETA-GAMMA  FILTER'S  SMOOTH  ACCELERATION  ESTIMATE. 

C 

A1 =0. 1 5 
A2=0.51 

B2=(A2**2)/ (2.0-A2) 

A3=0. 46864 
B3=0. 1 01 36 
6=0.000684 
C 

C  CALCULATE  THE  VARIANCE  THRESHOLDS  FOR  THE  ALPHA  FILTER  AND  THE  ALPHA- 
BETA  FILTER. 

C  V0UT1 VARIANCE  THRESHOLD  OF  THE  ALPHA  FILTER. 

C  V0UT2=VARIANCE  THRESHOLD  OF  THE  ALPHA-BETA  FILTER. 

C 

V0UT1  =  ( AT / 2 . 0-Al )*VRP 
V0UT1 =0. 6*V0UT1 

V0UT2= ( (2. *(A2**2)+B2*( 2. -3. *A2) )/ (A2*(4 . -B2*A2) ) )*VRP 
C 
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C  SHIFT  SMOOTH  POSITION  ESTIMATES  IN  SHIFT  REGISTERS 
C 

DO  5  1=1,  L2 

SP1 (L+l-I )  =  SPl(L-I) 

SP2{L+1 -I )  =  SPl(L-I) 

5  CONTINUE 

C  START  FILTERS  AFTER  TWO  MEASUREMENTS 
IF  (T.GF.0.2)  GO  TO  10 

C  INITIALIZE  THE  SMOOTHED  VALUES  OF  POSITION,  VELOCITY  AND  ACCELERATION 
SP1(1)=0 
$P2(1 )=0 
SP3=0 
SV1  =0 
SV2=0 
SA3=0 

C  COMPUTE  PREDICTED  POSITIONS  AND  VELOCITIES  FOR  THE  THREE  FILTERS 
10  PP1=SP1(1) 

PP2=SP2(1 )+TT*SV2 
PP3=SP3+TT*SV3 
PV1  =0 
PV2=SV2 

PV3=SP3+TT*SA3 

C  COMPUTE  THE  SMOOTHED  ESTIMATES  OF  POSITION,  VELOCITY  AND  ACCELERATION 
C  FOR  THE  THREE  FILTERS 

SP1(1)=(1.-A1)*SP1+A1*SP1(1) 

SP2(1 )=PP2+A2*(RADR-PP2) 

SP3=PP3+A3*(RADR-PP3) 

SV1=0 

SV2=PV2+(B2/TT)*(RADR-PP2) 

SV3=PV3+(B3/TT)*(RADR-PP3) 

SA3=SA3+C  G/(TT**2))*(RADR-PP3) 

CALCULATE  THE  MEAN  OF  THE  ALPHA  AND  ALPHA-BETA  FILTERS' 

SMOOTHED  POSITION. 

MSP1=0 
MSP2=0 
DO  20  1=1,  L 
MSP1  =MSPl-*-SPl  (I) 

MSP2=MSP2+SP2(I) 

CONTINUE 
MSP1=MSP1/L 
MSP2=MSP2/L 

CALCULATE  THE  VARIANCE  IN  SMOOTHED  POSITION  OF  THE  ALPHA 
AND  ALPHA-BETA  FILTERS. 

DO  30  1=1,  L 

D1 FI (I )=(SP1 (I)-MSPl )**2 
D1 F2( I )=(SP2(I )-MSP2)**2 
30  CONTINUE 
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VSP1=0 
VSP2=0 
DO  40  1=1,  L 
VSP1=VSP1+D1F1  (1) 

VSP2=VSP2+D1 F2( I ) 

40  CONTINUE 
VSP1 =VSP1/L 
VSP2=VSP2/L 
C 

C  IF  THE  VARIANCE  IN  SMOOTHED  POSITION  OF  THE  ALPHA  FILTER  IS  LESS  THAN 
THE  ALPHA  FILTER'S  VARIANCE  THRESHOLD,  THEN  OUTPUT  P0SITI0N=SP1 (1 ) 

AND  0  "PUT  VEL0CITY=0 
C 

C  IF  (VSP1.GT.V0UT1)  GO  TO  50 
SMP=SP1 (1 ) 

SMV=0 
RETURN 

IF  THE  VARIANCE  IN  SMOOTHED  POSITION  OF  THE  ALPHA-BETA  FILTER 
IS  LESS  THAN  THE  ALPHA-BETA  FILTER’S  VARIANCE  THRESHOLD,  THEN  OUTPUT 
P0SITI0N=SP2 ( 1 )  AND  OUTPUT  VELOCITY=SV2. 

IF  (VSP2.GT. V0UT2)  GO  TO  60 
SMP=SP2(1 ) 

SMV=SV2 
RETURN 

IF  THE  SMOOTH  POSITION  ESTIMATES  OF  THE  ALPHA  FILTER  AND  THE  ALPHA- 
BETA  FILTER  ARE  FOUND  TO  BE  INACCURATE  THEN  OUTPUT  P0SITI0N=SP3  AND 
OUTPUT  VEL0CITY=SV3. 

60  SMP=SP3 
SMV=SV3 
RETURN 
END 


END 


DATE 

FILMED 


DTIC 


